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ABSTRACT 
A model, called MAXPRESS, was constructed for the 
constrained optimization of the press cycle of the 


waferboard production process. Two techniques were used in 


the optimization routine. The first technique, the 
Hooke-Jeeves Direct Search Algorithm, allows the 
optimization of non-differentiable, non-linear, 


discontinuous, or undefined functions. The second technique 
employed was Everett's method of Lagrange multipliers. This 
method transforms a constrained optimization problem into an 
unconstrained optimization problem, and thus allows the use 
of the Hooke-Jeeves Direct Search Algorithm. 

MAXPRESS is a totally interactive program which prompts 
the user for all required input, and prints variable 
profit/shift and panel quality, as well as other parameters, 
at intermediate and final steps in the solution. Testing and 
evaluation indicated satisfactory model performance. Several 


recommendations are made for future research. 
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GLOSSARY 


feet per minute. 


furnish : The wood particles from which waferboard is 


IB 


MCF 


MOE 


MOR 


MSF 


O.D. 


psi 


formed. 

: Internal Bond; tensile strength perpendicular to the 
plane of the material; usually expressed in psi or Pa 
(pascal) (National Standard of Canada 1978). 

: Thousand (M) Cubic Feet. 

: Modulus Of Elasticity; the ratio of stress to strain, 
provided the stress does not exceed the elastic limit 
(and rupture the material); usually expressed in psi or 
Pa (pascal) (National Standard of Canada 1978). 

: Modulus Of Rupture; the maximum stress that a material 
can withstand before rupturing; uSually expressed in psi 
or Pa (pascal) (National Standard of Canada 1978). 

Thousand (M) Square Feet. 

: Oven Dry; term used in association with weight; O.D. 
weight is determined by drying the specimen at a 
temperature Of e052 C) meunty hima constant .Weignt is 
reached. 


pounds per square inch. 


waferboard : A type of board product made predominantly from 


large wood '‘wafers' (at least 30 mm. long), bonded 
together with a resin binder. Waferboard is a structural 
material and has been used extensively as a replacement 


for plywood. 
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3/8" equiv.: 'three eighths equivalent'; a base thickness 
for waferboard; most production figures are quoted in 


3/8" equiv. 


xi 


1. INTRODUCTION 

This dissertation describes the development of an 
optimization model for the waferboard press cycle. This 
model is the most important component of a larger project to 
develop an optimization model for the entire waferboard 
production process‘. The basic purpose for the development 
of a mill optimization model is to allow mill personnel to 
plan optimal management actions, especially in response to 
changes in the business environment, in a quick and 


inexpensive manner. 


1.1 Model Scope 

There appears to be three categories of questions which 
may be addressed with a waferboard production optimization 
model. The question(s) one wishes to address’ should 
determine the type of optimization model constructed. These 
questions deal with; 1) operating policy; 2) mill design; 


and 3) mill re-design. 


1.1.1 Operating Policy Questions 

The primary question affecting operating policy is how 
to achieve the most efficient operating policy with current 
mill equipment while meeting production and product quality 
objectives. Much of the equipment used in the manufacture of 
waferboard is highly specialized (and often custom built), 


‘Funding for this project was provided by the Forest 
Products Program of the Alberta Research Council. 
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and/or very expensive. Thus, the equipment set-up cannot be 
changed as readily as in some sawmilling situations (where 
there is typically a higher degree of equipment 
standardization and lower unit cost of equipment relative to 
waferboard manufacturing equipment). Therefore, a logical 
approach to the optimization of the waferboard production 
process over short time periods 1s to consider a current 
mill configuration as being fixed, and to establish the 
optimal operating policy for this configuration. 

Given these assumptions, it becomes evident that the 
efficient operation of the press is essential to an 
efficient mill operating policy. The press is the most 
important, and most expensive, piece of equipment in a 
waferboard mill. Operating policy must be centered on the 
press, and the pressing cycle, because this 1S where the 
production variables (eg. moisture content and resin content 
of the furnish, press time, etc.) ultimately have an effect 
Oneathe final product) in terms of both quality of product 
and’ tacesoraprodiictionss intthis contextea operating» opolicy® 
refers to the selection of a set of values for the 


controllable production variables. 


1.1.2 Mill Design Questions 

The most important issue in this category is how to 
design a waferboard production facility so as to achieve a 
given level of production with optimal efficiency. Mill 


design questions are long term in nature, and very little 
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can be considered as 'fixed", unlike the operating policy 
questions where the existing mill configuration is fixed. 
Efficiency implies a relationship between both costs. and 
benefits. Costs include the capital costs of the mill and 
related equipment, plus the capitalized expected costs of 
operating the mill over its life-span. Benefits include the 
Capitalized expected revenues gained from operating the mill 
over the same time period. This problem could probably be 
reduced to one of 'line balancing', perhaps using queuing 
theory, with consideration being given to likely future 
operating policies. The optimal mill design would change 
with different Operating policies. Thus, it iS important to 
consider flexibility as a requirement to the mill design 
question. This would involve the identification of key 
production variables which are likely to change in either 
Eves. Quantity, Or cost (or some combination of these), and 
the quantification of their potential impact on operating 
policy. 

Thus, it 1s clear that mill design questions are 
closely linked to operating policy questions. One must 
assume a mill design in order to establish an optimal 
eee policy. In the design of a mill, however, one must 
make assumptions about the likely (and, hopefully, optimal) 


operating policy which will be used. 
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1.1.3 Mill Re-Design Questions 

Mill re-design questions link the first two categories 
described above. Here, one is investigating the possibility 
of re-designing portions of an existing mill to achieve a 
more efficient operating policy and a better mill design. 
Also, one is addressing questions of equipment replacement 
necessitated by the failure of existing equipment. The 
capital costs of the alternate investment opportunities 
(i.e. between different pieces of equipment) are weighed 
against the capitalized profits (revenues less expenses) 


resulting from the re-design. 


1.2 Discussion 

The original mandate of this project was to 'develop an 
optimization model for waferboard production'. Initially, 
the author sought to develop an overall mill model which 
would allow one to answer all three types of question 
outlined above. A survey of available literature revealed 
very little work in the area of waferboard (or 
particleboard) mill optimization models. In fact, only two 
sources, Harpole (1979) and Balmasoff (1975), were found 
which deal with this topic at all. The computer programs 
presented in both papers are econometric simulation (not 
optimization) models. They do not consider the physical 
flows of materials within a mill, and can only provide very 


general answers to mill operating policy questions. 
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Numerous models exist for the Simulation or 
optimization of sawmills, or components of sawmills. Aune 
(1974) provides a good discussion of this type of mill 
model. Most such models utilize queuing theory to simulate 
the movement of logs and sawn lumber through the sawmill. 
These models determine the steady state conditions of the 
queuing system, and the user must resort to trial and error 
to test different mill configurations. This technique has 
some applicability for simulation of portions of the 
waferboard production process, where discrete individuals 
(such as logs or panels) move from station to station. It 
Still does not address the problem of optimizating either 
operating policy or mill design. 

A queuing model, coupled with an optimization routine, 
could be both feasible and practical for modelling some 
portions Are eee waferboard production process. Efforts were 
devoted to implementing an existing program (called DSMIN), 
developed by Carino and Bowyer (1979, 1981). DSMIN is a 
program designed for the optimization of some sawmill 
systems which can be represented by queuing models. The 
chief function of DSMIN is the optimization of maa design 
and re-design, with some application to the operating policy 
of sawmills. 

Soon after establishing DSMIN on the University of 
Alberta computer, it became apparent that it would be useful 
only for the simple queues of a waferboard production 


process. Queues which have batch arrivals (this is the case 
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when the press is loaded and unloaded), cannnot be 
accomodated with this program. More important, it became 
evident that neither DSMIN, or any other queuing model, 
could enable one to properly model the complex operating 
policy questions encountered in the waferboard production 
process. Unlike sawmilling systems, the waferboard process 
involves numerous, and complicated, management decisions 
about production variables which are related indirectly to 
throughput of individual pieces of equipment. These 
production variables include moisture content of the 
furnish, the amount of resSin applied to the furnish, panel 
density, wafer dimensions, press time, etc. (see Decision 
Variables - The Press Cycle in Chapter 3). 

It was then decided that the first stage in the 
construction of the overall waferboard production 
optimization model should be the development of a sub-model. 
The sub-model should adequately deal with the operating 
policy questions. Eventually it could be linked with a 
queuing model such as DSMIN (but able. to handle batch 
transfers as well), to anever all three categories of 
question about the entire mill. 

Operating policy questions must be adequately answered 
before one can properly deal with mill design and re-design. 
Operating policy greatly affects the throughput of mill 
equipment which subsequently has a large impact on _ the 
latter two categories of question. In addition, for existing 


mills, it seems that operating policy is more important than 
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mill re-design. Much of the equipment in a waferboard mill 
is very expensive (often custom built), and cannot, 
therefore, be easily replaced or modified to achieve optimal - 
throughput. 

It soon became apparent that the press cycle should be 
the focus for this stage of the waferboard production 
optimization model. It is in the press cycle that all of the 
production variables come into play to affect both the 
Guarelyeoreche=panelsproGuct "and the output vy ¥of*® -the**mrll 
(through their effects on press time). This dissertation 
describes the development of a model which optimizes the 
press cycle and, hence, the principle components of the 
operating policy of a waferboard mill. 

The next section of this dissertation contains a 
discussion intended to provide the reader with a basic 
understanding of the waferboard production process and the 
waferboard press cycle. After this, the two principle 
mathematical tools used to develop the model are discussed. 
These are the Hooke-Jeeves Direct Search Algorithm and 
Everett's method of Lagrange multipliers. This is followed 
by a complete discussion of the model, including a model 
user's guide, and an evaluation of model performance. A 
concluding chapter includes both general statements about 
the performance of the model and suggestions for additional 


research. 
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2. THE WAFERBOARD PRODUCTION PROCESS 

This discussion, and schematic diagram (APPENDIX I), 
are included EO provide the reader with a = basic 
understanding of the waferboard production process. The 
Weldwood (of Canada) Ltd. waferboard mill, situated near 
Slave Lake, Alberta, iS used aS an example of a typical 
waferboard mill. Elements of the production process at the 
Weldwood mill are described in the order of the production 
flow, from the woodyard, where green aspen logs are stored, 
through to the end product warehouse. The mill complex has 
been divided into three sections for this discussion; 
1) Harvesting and Woodyard; 2) Green Building, and; 3) Dry 
Building. Special reference shall be given to features of 
the production process which could have a direct bearing on 


production modelling efforts. 


2.1 Harvesting and Woodyard 

Weldwood hires logging contractors to harvest the aspen 
(Populus tremuloides) feedstock for the mill. Harvesting is 
conducted by the 'shortwood' method whereby whole trees are 
felled, skidded to a landing, and there bucked into 103" 
long logs. The logs are transported by truck to the woodyard 
adjacent to the mill, where they are unloaded. Wood in the 
Varo 1S) tOcated oneauarirst Ins erarst, Outi basis. 

The woodyard is a large reservoir of material for the 


mill. From a modelling perspective, this reservoir can be 
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viewed as an infinite supply of discrete individuals. 


2.2 Green Building 

The Weldwood mill is divided into two main buildings. 
The ‘green building’ (the wood has not been dried, hence, 
pomecia) nouses@ihat portion, of the yoroduction process /jirom 
log thawing and debarking, through to wafer drying. The ‘dry 
building' (the wood has been dried, hence, 'dry') houses 
that portion of the production process from dry storage to 
actual panel formation and warehousing. 

Logs are delivered in small bunches to _ the green 
building by a grapple loader. The logs are deposited into 
Gnew ore three loqg/ponds. These ilog) ponds) are constructed of 
reinforced concrete, ena Serve principally to thaw the wood 
in the winter, and to clean the wood during all seasons. 
Logs are forwarded down the length of the _ ponds by 
forwarding chains, and arrive at the green building entrance 
individually from each log pond. Logs from all three ponds 
may enter the boconGhee the same time, however, they are 
transferred to a Single conveyor and are forwarded to the 
debarker as a single queue. 

Logs are then fed through a rotating ring debarker, one 
at a time, where most bark is removed. Feed rate through the 
debarker can be varied, depending on the condition of the 
logs and bark (eg. how large the logs are, degree of 


thawing, etc.). 
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The next stage iS the slasher deck. Logs are _ fed 
broadside, one at atime, into this set of circular blades 
ancdeare, cutwintowsrour 425") "bolts, sThese ss bolts are’ then 
forwarded to the waferizers. 

Three disc waferizers are used at the Weldwood mill. 
Normal operation has two cae Ups mone) junit” down for 
maintenance (blade change, etc.). The waferizers accept 
bolts from the slasher deck, and feed them broadside into 
the knives on the rotating discs. Each waferizer has its own 
input queue of bolts which is fed from the single queue 
coming out of the slasher deck. Feed rate into the 
waferizers varies depending on the desired wafer thickness; 
the thicker the wafer, the faster the feed rate. Output from 
the waferizers is in the form of a flow of wafers. 

The flow of wafers then passes through one of two 
rotating drum type 'green' screens, which separate the green 
'fines' (small particles unsuitable for the manufacture of 
waferboard) from the green wafers. 

The flow of green wafers then passes to one of three 
green storage bins which feed the wafer dryers. These bins 
provide a buffer for the production process. From the’ green 
storage bins, the wafers flow continuously into three 
natural gas-fired dryers, and then flow out of the_ green 
building. The time wafers spend in the dryers depends on the 
moisture content of the green wafers, the desired moisture 
content of the dry wafers, and the temperature in the 


dryers. The dryers at the Weldwood mill are each rated at 
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7000 O.D. pounds of wafers/hour. 


2.3 Dry Building 

The flow of dry wafers passes out of the green building 
taco eehe. dryssvorage, bini-of .the ) dry obuidding.: ‘The. dry 
Storage bin provides another buffer to the production 
process. Wafers pass through the dry storage bin on 
approximately a First In - First Out basis. After leaving 
the dry storage bin, the flow of wafers is separated, on the 
basis of wafer size, by a Screen and vacuum pick-up, into 
two different flows. Each flow then passes through one of 
two rotating drum type dry screens (same function as the 
green screens). Large wafers (for the faces of the 
waferboard panels; comprising approximately 40% of the 
total) are fed through one screen, while small wafers (for 
the core of the waferboard panels; comprising approximately 
60% of the total) are fed through the other. 

Following screening, the two types of wafer pass into 
two surge bins, one bin for each type of wafer. Each surge 
bin mixes the flow of wafers passing through it, in order to 
achieve uniform bulk density. The two flows then pass _ into 
two rotary drum blenders where wax and phenolic resin are 
applied to the face and core wafers. The wafer flows’ then 
enter the press line. 

The first components in the press line are the forming 


heads. Wafers flow into four forming heads; two for the face 


Ue 


of the panel, and two for the core. Cauls (steel sheets; 
approximately 16'x4'), overlapping at their ends, pass under 
the forming heads where wafers are dropped to forma 
continuous mat with random wafer alignment. Caul speed, and 
the rate of wafer flow onto the cauls, vary with different 
mat thicknesses for the production of different panel 
thicknesses. Cauls exit from the forming heads and are 
separated from each other. 

A new queue is formed at the next stage, the 24-opening 
press loader, where cauls, with their wafer mats, are 
individually entered until all 24 positions are filled. The 
entire batch of cauls is =hen simultaneously transferred to 
the 24-opening press. The wafer mats are then pressed under 
high temperature and pressure (approximately 400°F. and 800 
psi). Press time (usually ranging from 3 to 7 minutes) can 
vary depending on a number of parameters including panel 
quality desired, panel thickness, moisture content of the 
mat, etc. All 24 cauls are then unloaded simultaneously from 
the press to the 24-opening press unloader. From there, each 
caul is unloaded individually and separated from the pressed 
defers mat. The pressed mats are rotated through a board 
cooler, and then pass through two sets of trim saws. The 
resulting 4'x8' panels are graded, bundled, and stored or 


shipped. 
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(2.4 Discussion 

The overall waferboard production process can be 
represented by a system of queues and flows (APPENDIX I). 
Queues consisting of logs exist at the debarker (single 
server) and the slasher deck (single server) while several 
(one, two or three) queues, made up of bolts, occur at the 
waferizers (multiple servers). A flow of wafers occurs from 
the waferizers to the forming heads. At the forming heads, 
another queue is formed, the elements of which are cauls. 
The cauls arrive at the press loader, and leave the press 
unloader individually, but are loaded and unloaded from the 
press in batches of 24. After being separated from the 
cauls, the pressed mats form queues at the panel cooler and 
at the trim saws. Grading and sorting could be viewed as 
part of the trim Saw operation or aS a separate queue. 


Panels are moved to the warehouse in batches. 
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3. THE MATHEMATICAL MODEL 


3.1 Assumptions 
The following basic assumptions have been made when 
formulating the mathematical model: 

1. The overall goal for the waferboard mill is to maximize 
profits. 

2. The mill must meet, or exceed, a single minimum level of 
Guality for its waferboard production. In other words, 
there iS no option for the production of a variety of 
panel grades which can then be sold for different 
prices. However, a variation of the standard application 
of the model can accomodate the case of varying panel 
qualities (refer to Economic Significance of the 
Lagrange Multiplier in Chapter 5). 

3. All of the waferboard produced which meets the minimum 
quality level can be sold at a specified factory gate 


price (this price is provided by the model user). 


3.2 Decision Variables - The Press Cycle 

Numerous variables in the waferboard production process 
affect final panel properties and mill revenues and costs 
directly, and through interactions with each other. These 
variables include resin amount, type and distribution among 
layers of the wafer mat, moisture level and distribution in 


the wafer mat, panel density, press time, temperature, 
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pressure and ltrate:-iof, closure; particle geometry and 
Orientation (if any), wood species, special additives, etc. 
(Maloney 1977). Many of these variables may be controlled, 
and could be considered decision variables for an 
Optimization of the waferboard production process. In 
addition, these variables ultimately have their effect on 
both board quality, and mill production rate, when the wafer 
mat iS pressed. Therefore, it becomes apparent that the 
press cycle is the most important part of the waferboard 
production process. Not only are final panel properties 
determined here, but in addition, press time will determine 
mill production rate and hence, revenue. 

Four controllable production variables were selected as 
decision variables for inclusion into the mathematical model 
developed for this study. Operating policy is defined by a 
set of values for these four variables; 

1. resin content of the wafer mat, 

2. press time, 

3. moisture content of the wafer mat, 
4, nominal panel density. 

These variables were selected because of their 
importance to panel quality and mill profits, and because 
data concerning their effects was available in the 
literature (refer to Constraint Data Variable in Chapter 6). 
It is important to note that other variables could, and 
probably should, be included into this model if sufficient 


data were obtained. The list of variables presented above is 
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by no means exhauStive, nor are the variables included 
necesSarily the most significant ones in the waferboard 
product i.onweprocessssTheyatdo, however, tillustrateb«model 


functioning very well. 


3.2.1 Resin Content 

As resin content in the wafer mat increases, all 
Strength properties increase (Kelly 1977, Maloney 1977). In 
general, however, these strength properties increase at a 
decreasing rate. 

Resin accounts for a major part of total mill expenses. 
The selection of potsned resin content, therefore, is very 
mmpontant tro miklaprofptability,;sA mill optimization’! tmodel 
should be constructed to allow the selection of the optimal 
resin content in response ele) resin price changes. 
Furthermore, most resins currently in use are synthetic in 
nature (produced from petroleum). It appears that resin 
prices will be closely associated with oil prices and could 
be subject to significant price increases in the future, as 
they have in the past (Maloney 1977). 

Many different types of resin are currently available. 
Each has its aan physical properties (for example, curing 
time) and cost. In addition, as synthetic resin costs rise, 
nore work 1s devoted to the development of new resins which 
can be manufactured from natural products (Maloney 1977, 
Dolenko and Shields 1980). Because each resin has its own 


properties, a reSin which is cheaper to purchase may not 
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necessarily be cost efficient. Such a resin could require 
longer press times for proper adhesion. This might result in 
a greater cost (through lost production) than the benefit 
gained through the lower resin prices. It is important to be 
able to properly evaluate new, or different, resins with a 


mill optimization model. 


3.2.2 Press Time 

Press time affects the rate of mill production (hence, 
mill profits) and final panel quality. In general, the core 
temperature of the wafer mat must reach the level required 
to cause the resin to cure, without subjecting the panel 
surface to excessive temperatures which result in 
degradation (Kelly 1977). 

In the actual operation of a waferboard mill, there 
appears to be ae relatively direct trade-off potential 
between resin content and press time (James 1981). Within 
limits, press time can be reduced if resin content is 
increased. Both panel quality and mill profits will be 
greatly affected by this trade-off. This is an important 
operating policy decision which is faced continually by mill 
managers. A mill optimization model should be capable of 


addressing this problem. 
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3.2.3 Moisture Content 

Moisture content of the wafer mat directly affects heat 
teansferetonthet@core jand, 6therefore vwitsetatecriticalsrfactor 
for both press time (therefore, mill profits, as above) and 
panel quality (Kelly 1977, Maloney 1977). Surface moisture 
evaporates and travels to the core when the mat is pressed, 
thus allowing quicker heat transfer than would otherwise be 
possible (Kelly 1977). However, excessive moisture requires 
increased press time to allow for adequate evaporation 
through the panel edges, and to compensate for any adverse 
effects the moisture may have on resin curing (Kelly 1977). 

Moisture decreases the compressive strength of wood and 
helps create a high density gradient in the finished panel 
product; the core tends to have lower density compared to 
the panel faces (Kelly 1977, Maloney 1977). This results in 
high bending ae ee parallel to the board surface (eg. 
modulus of rupture - MOR, and modulus of elasticity - MOE), 
and low tensile strength perpendicular to the plane of the 
panel (eg. internal bond - IB). 

The level of moisture content in the furnish is 
directly related to fuel consumption in the gas eres 


which is another component of total mill expenses. 
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3.2.4 Nominal Panel Density 

Panel density depends on the amount and species of wood 
used, resin, moisture, and additives content, and pressing 
methods (Kelly 1977, Maloney 1977). In general, strength 
properties will increase as panel density increases, 
provided press time is adequate (Kelly 1977, Udvardy 1979). 

Nominal panel density was included as a decision 
variable for this project because the Weldwood mill at Slave 
Lake uses 'Sstops' in its pressing routine. With the use of 
stops, variation of panel density is primarily a function of 
the mass of furnish in the wafer mat. This variable might be 
re-defined better as the mass of furnish on each caul 
entenanggnthepapressye In gany*icasey. thiswavariable helps 
illustrate the use of the optimization model. 

Panel density also has a direct effect on the costs of 
production. As more furnish is added to the wafer mat to 
increase panel density, wood costs/panel increase, as do 


resin costs, fuel costs, etc. 


3.3 The Optimization Model 

With the removal of marketing considerations, as 
assumed in Section 3.1, the basic problem addressed by this 
study is the optimization of a payoff subject to a resource 
constraint. In particular, for the waferboard press cycle, 
the problem has been defined as the maximization of profit 


while meeting a minimum level of panel quality: 
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MAXIMIZE: 


PROFIT = REVENUES - EXPENSES 
SUBJECT TO: 


QUALITY 2 MINIMUM QUALITY 


WHERE: 

PROFIT = £,{Resin Content, Moisture Content, 
Press Time, Panel Density} 

QUALITY = f,{Resin Content, Moisture Content, 


Press Time, Panel Density} 

In the model, profit is measured in dollars/8-hour 
shift. Panel quality 1S measured as internal bond in psi. 
Other measurements of panel quality could be used (eg. MOR 
or MOE) instead of, or in addition to, internal bond (refer 
to CONCLUSION AND RECOMMENDATIONS in Chapter 9). 


Lff assumptions, i@suchy tasvelinearitypuscontinulty; oor 


differentiability, could be made about the objective 
MunGuLon mo )juevana. *theseconstraint. “fLunctiom (Ci. )4 then 
classical optimization techniques such as linear 


programming, or differential calculus, could be employed to 
solve the constrained optimization problem. However, these 
assumptions cannot be made about f, or f, since almost 
nothing can be guaranteed about the nature of either 
relationship. In fact, the two relationships may not even be 
defined as functions (see Design Considerations in Chapter 
6). The methods employed to solve the constrained 
optimization problem must therefore be capable of accepting 


the objective and constraint relationships in either 
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EUNnctIOnal .form@ion-astimatrices of “datas 

These characteristics led to the selection of two 
relatively non-classical optimization techniques, the 
Hooke-Jeeves Direct Search Algorithm, and Everett's method 
of Lagrange multipliers, to accomplish the constrained 
optimization. Everett's method of Lagrange multipliers 
transforms the original, constrained optimization problem 
into an unconstrained problem. The optimal solution of the 
unconstrained Lagrange problem can then be found using the 
Hooke-Jeeves Direct Search Algorithm, without assumptions 
about the nature of the functional relationship. The 
Hooke-Jeeves Direct Search Algorithm is described in Chapter 
4 of this dissertation. Chapter 5 contains a discussion of 
Everett's Method of Lagrange Multipliers. Details of the 
computer version of the model are then discussed in Chapter 
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4. HOOKE-JEEVES DIRECT SEARCH ALGORITHM 

The Hooke-Jeeves Direct Search Algorithm is a routine 
(called the 'pattern search routine') for optimizing 
(minimizing or maximizing) an unconstrained function S(¢) of 
Snieyre Monel var vables) p= tp 4h 0555'055 Cone? go ) (Hooke ‘and 
Jeeves 1961). | 

The values of ¢ can be interpreted as points in 
N-dimensional space. A 'move' may be thought of as a vector 
projection in this space and is defined as the procedure of 
going from a given point (representing a particular value of 
S(g)) to another point in the same N-dimensSional space. A 
move is termed a 'success' if the value of S(g) is improved 
(i.e. if S(g) decreases in a minimization problem or 
increases in a maximization problem). A move is termed a 
"failure' otherwise (Hooke and Jeeves 1961). The pattern 
Search routine makes two types of moves; 1) exploratory 


moves; and 2) pattern moves. 


4.1 Exploratory Moves 

Exploratory moves acquire information about’ the 
behaviour of the function S(g) solely by their success or 
failure (Hooke and Jeeves 1961). This information is 
utilized to establish a probable direction for a_ successful 
pattern move. Exploratory moves are achieved by changing the 
coordinate values (¢), one at a time, and comparing the new 


value of S(g) to its previous value. Each coordinate (9g) is 
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first raised by some arbitrarily chosen step size, and S(g) 
Vs) evaluated. §If (the move is’.a failure, ‘the original 
coordinate (g) is lowered by the same step size and S(¢) is 
re-evaluated. Each time a successful move is achieved, the 
Original set of coordinate values is re-set to include the 
coordinate value (¢g) by which the eer eral move was 
achieved. The routine then moves on to the next coordinate. 
This procedure 1s carried out sequentially for all 


g. (i=1,2,3 ...N) (Figure 1) (Hooke and Jeeves 1961). 


4,2 Pattern Moves 

Pattern moves are designed to utilize the information 
acquired by the exploratory moves about the behaviour of the 
function S(¢). A base point is defined as a point 
(represented by some particular set of coordinate values 
% (i=1,2,3 ...N)) from which a pattern move is made. Thus, 
pattern moves may be viewed as proceeding from base point to 
base point, with the pattern move from a given base point 
duplicating the combined moves from the previous base point 
(Hooke and Jeeves 1961). This is accomplished by changing 
all coordinates (9, (i=1,2,3 Ne uN) DY any seamount. Vequalu)co 
the difference between the present base point, and the 
previous basepoint (Hooke and Jeeves 1961). As ae result, 
once a pattern is established, the size of the pattern moves 


will increase, resulting in an acceleration in that 


direction (Phillips et al. 1976, Hooke and Jeeves 1961). 
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Flow Diagram for Exploratory Moves of Pattern Search Routine 


Ax 0G 


Legend 


The variables ¥, 6 and © are points in a 
K-dimensional space; the rest of the variables are 
unidimensional. 


the previous base point 

the current base point 

the base point resulting from the current move 
- the test base point 


s(¥) the functional value at the base point 
s(@) the functional value for this move 


the functional value before this move (usually, the 
smallest value so far attained by the set of 
exploratory moves) 

current step size 

“minimum” step size 

reduction factor for step size, p <1 

one of the coordinate values for @, k = 1,2,...,K 
number of coordinates for the points 
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The argument for pattern moves iS intuitive in nature; if a 
set of moves was successful in the past, it is likely to be 
successful again in the future (Hooke and Jeeves 1961). The 
pattern move procedure is clarified by the following 


discussion of the overall pattern search routine. 


4.3 Pattern Search Routine 
A detailed flow diagram of the pattern earch routine 
is) shown, ingkigure2 and tsegraphical analogy’ is shown in 
Fngqure 3.. Using) function (minimization. form illustrative 
purposes, the following steps are performed in the pattern 
search routine (Hooke and Jeeves 1961): 
Sipe) by, 
- Start at the ‘current base point' (arbitrarily chosen 
for initial iteration): 
- Make exploratory moves. The base point resulting from 
these moves is the 'test base point’. 
- Is the functional value at the test base point below 
that at.thescurnent base point? 
TehoemGOrcO sot ae 2. 
SNOT GOawCO To hE rads. 
STEP 23 
- Re-set base points; 
‘current base point' becomes 'previous base point'; 
"test base point' becomes ‘current base point’. 


- Make pattern move by doubling the coordinates of the 


* ‘ . 
Tt yretsa ») 
@ tt Vaver ba 
om od qienit «! whet ood. 1 


ia Vee 


p&iwot fas ete yao ae ed, 
- 


onli Ar ae 


i260 oerns wics2aq edt Fo expel wot 
n? neota a2 voéigta lesldgase @ ot V6 i. Pb rs 
t seen yi! [) po. tettetde acy ae 7 

iy eT Lee’ Sg S46 aeate paewe that aio 
(CCE! Revet) tae gamit) ila 


a 


ia. 
send net Tap gt3 18: 


a 
SWotapinet pei 
igluea oAled bes 2fY .2svnil> eres ee aoe < 
i . ’ 
, Oftag. ees *Osl poh ei 2590R ceil 
if =e in iog Lp, 2Ja99 Sent If Bece iatot tong 3a3 at - . 
C3ntuq eset tosivys edly ae eds’ 
2 GATS of OD Taa¥ =. 
-~ Gitc o2 68 5 O8 =) 


a» 
i, 
* 
c 
—_ 
mi 
, 
~~ 
a 
- 
= 
~ 
™ 


 teiniog seed 2 
intog cas cvolverg? sameted.gniog ser jowrigs" 
‘cniog sned sasiahe’ sence “1nheq sted se 


ei 3m) gatani 5360s. ats 
: ra 


Figure 2. 


+(2¢) - 6 
s=-s() 


Flow Diagram for Pattern Search Routine 
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Legend 


The variables ¥, 6 and © are points ir a 
K-dimensional space; the rest of the variables are 
unidimensional. 


the previous base point 

the current base point 

the base point resulting from the current move 

- the test base point 

the functional value at the base point 

the functional value for this move 

the functional value before this move (usually, the 
smallest value so far attained by the set of 
exploratory moves) 

current step size 

“minimum” step size 

reduction factor for step size, 9 <1 

one of the coordinate values for $, k = 1,2,...,K 
number of coordinates for the points 


exploratory search routine 
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test base point (now the same as the current base 

point), and subtracting the coordinates of the previous 

base point. This is the test base point. 

- Make exploratory moves. The coordinates of the test 

base point may, or may not, be changed by these moves. 

- Is the functional value at the test base point below 

that at the current base point? 

Bo oe GOmCOmorT bri .. 

SRN ws GO 1bOe ore 1%, 
STEPS « 

Zar seeScTeDp, «size» for exploratory moves small enough 

(compared to an arbitrarily chosen minimum) ? 

aN Sc) ho LOOP. 

- NO : Decrease step size; Go to STEP 1. 

Hooke and Jeeves (1961) suggest that all pattern moves 
be immediately followed by exploratory moves (STEP 2), 
before testing the pattern move for success. Their rationale 
is as follows: 
Because more progress towards optimizing the function S(@) 
is made with pattern moves than with exploratory moves, it 
is desirable to retain pattern moves where possible. A 
pattern move which otherwise would have failed can sometimes 
succeed and, hence, be retained, if exploratory moves are 
made after the pattern move and this esis tested for 
Success. Thus, the success or failure of a pattern move 
becomes irrelevant to the mechanics of the routine because 


in either case, exploratory moves are made immediately 
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TOUTOyiungethespaLrern i moveg Thi sestrategyersemostilikely, | to 
Succeed when the pattern is first being established (Hooke 


and Jeeves 1961). 


4.4 Discussion 

There are several practical considerations to the 
Hooke-Jeeves Direct Search Algorithm which affect potential 
applications and are, phonercres discussed here. 

A primary advantage to using the Hooke-Jeeves Direct 
Search Algorithm is that the objective function to be 
optimized does not. have to be regular, continuous, or 
differentiable, nor does it have to be explicitly defined 
(Carino and Bowyer 1979, Phillips et al. 1976, Hooke and 
Jeeves 1961). This is a particularly useful feature for 
applications such as the waferboard production model 
presented in this thesis. In this case, the objective 
function is a combination of an explicitly defined profit 
function and a constraint relationship represented by 
Br Screve data in la Sere: Formate. adareLon., the 
algorithm lends itself well to use on computers, since it 
uses repeated arithmetic operations with simple logic (Hooke 
and Jeeves 1961). 

A disadvantage to the Hooke-Jeeves Direct Search 
Algorithm is that the global optimum will not always be 
found. With some objective functions (ie. those defining a 


non-convex feasible region) it is possible, even probable, 
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that some of the solutions produced by this algorithm will 
be local optima only. In such cases, it appears that the 
only practical solution is to start the search from several 
arbitrarily chosen “points “and compare the resulting 
solutions. This will not, however, guarantee finding the 
global optimum. It only increases the chances of such an 
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5. EVERETT'S METHOD OF LAGRANGE MULTIPLIERS 

Everett's method of Lagrange multipliers is useful for 
optimization subject to constraints, especially in problems 
where discontinuous or non-differentiable functions must be 
optimized (Everett 1963). This property of Everett's method 
is particularly useful to some applications, such as the 
waferboard production model described in this dissertation, 
where little is known about the nature of either the 
function to be optimized or the constraints which limit 
possible solutions (eg. functional form, linearity, etc.). 
Clearly more conventional techniques for constrained 
optimization, such as linear programming, could not be used 
“3 BES Situation. 

The following discussion of Everett's method of 
Lagrange multipliers will deal primarily with the practical 
applications and problems associated with this technique. 
However, some discussion of the theoretical justification 
for Everett's method is warranted, and follows in the next 


section. 


5.1 Main Theorem 

Using the terminology of Everett (1963) the main 
theorem shall be discussed in terms of the optimal 
allocation of limited resources. In other words, the problem 
is the maximization of a payoff function subject to given 


constraints. The waferboard production process provides a 
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good example. In this process the objective is to maximize 
pootius yn Subjectiegco s¢qualityyconstraints, onetheRwaferboard 
panel product. Quality may be measured by internal _ bond 
(IB), modulus of rupture (MOR), modulus of elasticity (MOE), 


SEC. 


5.1.1 Definitions 

The following definitions are used in the discussion of 
the main theorem, and are derived from Everett (1963) with 
one major change; only one resource function and constraint 
is considered?: 

1. x: The decision variables (i.e. resin content, press 
time, moisture content, panel density). 

2. S: The set of possible strategies. (eg. the possible 
combinations of the decision variables to produce 
waferboard of various quality levels, at various 
production rates). 

3. H: A real valued function called the 'payoff function’ . 

.4, H(x): The payoff which occurs as a result of employing 
the strategy os (eg. could be expressed as 
Purohit /S-houmeshitt’) 

570% tA real svalued function jcalled Be "resource function’. 

6. C(x): The resource expenditure (panel quality is the 
resource considered in the waferboard press cycle) 


required to gain payoff; occurs as a result of employing 


2 The waferboard production model allows for only one 
constraint (internal bond); see Constraint Data Variable. 
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the strategy x¢S. In other words, some resource must be 
given up in order to gain payoff. This could be thought 
Of asimean mm .expenditure ition! “loss” pict quality: leg. 
measured by internal bond) in the waferboard panels in 
order to achieve more profit. 


7. c: The maximum resource expenditure allowed. 


5.1.2 Discussion 
'The problem then becomes: 
MAXIMIZE: 
H(x) over all (x¢S) 
SUBJECT TO: 
CaS 
The main theorem follows (for a proof of the main 
theorem, see Everett (1963)): 
1. A is a nonnegative real number, 
2. x*¢€S maximizes the function, 
HUx)-AC (x) over all x¢S, 
Grex ema xim: zesuH Ux) mover all, choses x,s) suUChsehatsG SVC x1). 
In other words, for any enouce of nonnegative 4, if an 
neoneer ened maximum of the Lagrange function (item 2, 
above), can be found (x* is the strategy which produces this 
maximum), then this strategy will also produce the maximum 
of the constrained function with constraints equal to _ the 
amount of the resource expended in achieving the 


unconstrained solution (Everett 1963). 
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Using the main theorem, one can arbitrarily choose a 
non-negative A and find the maximum of the unconstrained, 
Bagrange function. — This ’'Solution will” alsomiprovide the 
maximum of the Omlginale, = Fconstuained. func oon, eew1 oh 
constraints equal to the amount of resources expended in the 
unconstrained solution. It iS important to realise, however, 
that the choice of A is completely arbitrary, and different 
choices of A will generally lead to different resource 
levels. If one is interested in maximizing a payoff function 
(such as gross revenue/8-hour shift) while not exceeding a 
given resource expenditure (eg. loss of internal bond below 
some minimum allowable level), it is likely that the value 
of t will have to be adjusted by trial and error until the 


desired constraint level is achieved (Everett 1963). 


5.2 Economic Significance of the Lagrange Multiplier 

The Lagrange multiplier (2) also provides some 
information about the cost of the constraint. By the Lambda 
Theorem (for a discussion and proof of this theorem, refer 
to Everett (1963)) it can be shown that at any optimal 
solution, A represents the marginal value per unit of 
resource. In the waferboard production case, 2% would 
represent the profit gained by the ‘expenditure’ of one 
additional unit of quality (eg. loss of one unit of internal 
bond). This feature could have implications for a mill's 


operating policy because as resource costs change, so does 
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the cost of achieving a given level of panel quality. In 
acdition,, if, Oonewihas “the option “of. marketing. several 
Giiperent grades OL product, /iticould’ prove, optimal to vary 
the quality of the product in response to changes in 
resource prices and selling prices for the panel product. In 
other words, if the situation arose where a waferboard mill 
could actually make a higher profit by producing different 
quality panels, Everett's method of Lagrange multipliers 
could help to identify the optimal quality level of panel to 


produce. 


5.3 Gaps or Inaccessible Regions 

The use of Lagrange multipliers does not ensure that a 
solution will necessarily be found he all problems, 
however, a solution found by their use is guaranteed to be a 
true solution (Everett 1963). There is no guarantee, 
however, that some constraint levels will not be generated 
by any values of A. These are are termed inaccessible 
regions (also called 'gaps') (Everett 1963). 

A gap can be recognized by abrupt dicontinuities in the 
resource levels generated as & is continuously varied (i.e. 
for two levels of 4&4 which are very close, the corresponding 
resource levels generated are considerably different). 

The basic cause of gaps iS nonconcavity in the function 
of optimum payoff versus resources expended (Figure $4) 


(Everett 1963). The Lagrange method will succeed in 
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Producing rallysolutionss mnvieche: concave seregrons yeof Aithis 
hunct toni, and will fail in all non-concave_ regions. 
Solutions in the inaccessible regions cannot be obtained by 
a simple application of the Lagrange multiplier method, and 


must be sought using different techniques (Everett 1963). 


5.4 Method for Handling Gaps 

Several methods for handling gaps are presented by 
Everett (1963), but only one of these methods, the most 
applicable to the waferboard production model, will be 
discussed here. 

Consider the standard case where one is seeking to 
maximize some payoff function subject to constraints. If the 
decision variables are discrete (indeed, as will be 
explained later, they are discrete in the waferboard 
production model), then nearly optimal solutions can be 
produced by deliberately deviating slightly from the optimal 
combination of decision variables which last maximized the 
Lagrange function. This can be achieved by examining all 
possible combinations of the decision variables + one step 
(or more) away from’the last optimal combination (C) of the 
decision variables (i.e. the last combination which produced 
an optimum solution using the direct application of 
Everett's method). The Lagrange solutions are generated for 
all of these combinations, and are then subtracted from the 


Lagrange solution obtained with C. The resulting deviations, 
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with their corresponding combinations of decision variables, 
are then ordered from lowest to highest. A member of this 
ordering is 'dominated' if a preceding member provides more 
payoff for equal, or less, resource expenditure. Dominated 
members can then be dropped out of the ordering, since they 
are not optimal solutions. Remaining members of the ordering 
define acceptable strategies for eheeoonpe ation of decision 
variables. These are represented by corner points ('OTHER 
SOLUTIONS') in the non-concave regions of the optimum payoff 
vs. resources expended function of Figure 4. If a minimum 
constraint level must be achieved (such is the case in the 
waferboard press cycle, where minimum panel quality 
Standards must be met), the undominated combination of 
decision variables which provides a constraint level closest 
to, and greater than, the minimum constraint level allowed 
can be selected. This will be the best possible solution 
because no other members of the ordering will provide the 
same, or more, payoff for equal, or less, resource 
expenditure. In other words, for the waferboard press.cycle, 
the undominated member (combination of decision variables) 
of the ordering aon provides a level of panel quality 
closest to, and greater than, the minimum specified level of 
panel quality, is selected. No other member of the ordering 
will provide the same, or more, panel quality (equivalent to 
the same, or less, expenditure of panel quality) with the 


Same, Or more, profit. 
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In some cases, such as the waferboard press cycle, 
where little is known about either the payoff or constraint 
relationships, one could guarantee that the best solution 
(i.e. undominated member of the ordering closest to the 
minimum quality level) is found only if all feasible 
combinations of the decision variables are examined. 
COOVOus Ly sethiismi SenOt- practical ..ueime fact, sstinius tstrategy 
would make~ the use of Everett's method of Lagrange 
multipliers obsolete, because all possible solutions would 
be generated directly. In this case, justification for 
examining combinations of decision variables within a 
limited range of C is intuitive in nature, and is supported 
by empirical results. It seems logical, in most cases, to 
assume that solutions (i.e. combination of decision 
variables) occurring in gap regions will be close to _ the 
last known solution which bounds the gap. This assumption 
was tested, and found to be valid, for gaps produced by the 
objective and constraint relationships developed for this 
Study. All possible combinations of decision variables were 
generated, and the best combination was selected for the 
resource level desired. For all cases tested, this result 
was the same as the result obtained using gap search 
techniques directly. 

The same justification could be presented for the use 
of gap search techniques around an optimal solution obtained 
directly with Everett's method, but which is not necessarily 


at the edge of a gap. If a resource level generated with a 
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particular value of A is close to the desired resource 
level, a direct application of gap search techniques could 
provide the optimal solution quicker, and easier, than 
further manipulation of AR Once again, the only 
justification for this procedure is intuitive in nature, and 
Supported by empirical results. Since the underlying aim of 
this» project, was to construct a ‘useful; practical tool,  1t 
was felt that these gap search procedures should be 


available in the waferboard production model. 
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6. THE COMPUTER MODEL - MAXPRESS 
The model developed for this dissertation has _ been 
named MAXPRESS?, an opiehak pen fee model for the press and 
pressing cycle of a typical waferboard production facility. 
MAXPRESS maximizes variable profit/8-hour shift while 
meeting minimum panel quality specifications (measured as 


internal bond in psi). 


6.1 Design Considerations 
MAXPRESS has been constructed with the following 
conditions and desired features in mind: 

1. MAXPRESS should be able to accept non-differentiable, 
non-linear, and discontinuous objective (payoff) 
functions. The author felt that a payoff function could 
be defined, but that very little could be guaranteed 
about its properties. In addition, MAXPRESS should be 
constructed so that it could be easily modified to allow 
for the situation where the objective function could not 
be defined, and where the payoff information was 
provided in the form of discrete data for various 
combinations of decision variables. 

2. MAXPRESS should accept constraint information in the form 
of discrete data for various combinations of decision 


variables. This condition was recommended by Dr. Lars 


3 MAXimization of the PRESS cycle. 
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Bach* and was verified by an examination of waferboard 
quality relationships in the literature. Good functional 
relationships between waferboard quality (eg. internal 
bond) and production variables (eg. resin content, press 
time etc.) appear to be non-existent in the literature. 
This might be due to the complexity of the relationships 
involved. It does appear, however, that fairly good data 
can be obtained from the literature (certainly, any 
forest products company wishing to use a model such as 
MAXPRESS could provide this data) in the form of 
repeated measurements (or averages of such measurements) 
of quality for various combinations of decision 
variables. A major problem with such data obtained from 
the literature 1s the wide variety and incomplete 
description of testing methods used. 

Furthermore, anomalies could occur in a particular 
waferboard mill which might Significantly alter 
quality/decision variables relationships. In such a 
case, it 1s desirable to have a model which could use 
Giatitys data vafnomestrial production runs, or from an 
"educated guess' of mill personnel, without trying to 
develop predictive functions using techniques such as 
linear, or non-linear, regression. 

3. MAXPRESS must be ng optimization model, not a simulation 


model. The ultimate goal, and mandate, of this project 


* Program Manager, Forest Products Program, Alberta Research 
Council and member of the author's examination committee. 
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LSMCOMIPLOdUCem@ an wrODtLimi Zationm model pefor= tan. sentirre 
waferboard production facility which can be used by mill 
personnel. 

4. MAXPRESS should be easy to use. For example, waferboard 
production personnel at a mill (eg. the quality control 
manager) with no formal training in computer programming 
Or operations research should be able to use the model 
to easily answer particular questions about the mill's 


operating policy. 


6.2 Model Construction 

MAXPRESS, the model constructed, meets alle “of the 
requirements described above. It iS an optimization model 
designed to accept any type of explicitly defined objective 
(payoff) function, and quality relationships in the form of 
discrete,. tabular data. Optimization is achieved by using 
the Hooke-Jeeves Direct Search Algorithm (pattern search 
routine) in combination with Everett's method of Lagrange 
multipliers. The pattern search routine was modified to 
allow discrete constraint inputs to the model. As a result, 
the step size for exploratory moves must equal the interval 
size of the constraint matrix, and no reduction in step size 
is allowed. No limiting assumptions are required for either 
the objective (payoff) function or the constraint (quality) 


data variable. 
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MAXPRESS is a totally interactive computer program 
which Promptssp nes users, for Vall arequired inputs; and 
automatically provides the required output. The model user 
does not require any special knowledge about computers or 
computer programming, other than how to access his_- own 
system. Some understanding of the role of the Lagrange 
multiplier (called LAMBDA in the model) is required (see The 
Role of LAMBDA in Chapter 7). 

MAXPRESS was formulated in APL (A Programming Language, 
Gilman and Rose 1976). All major routines are composed of 
different subordinate functions which are called by a main 
function as required. This modular composition, 
areeAatee ge of programs written in APL, permits easy 
substitution of various functions or the addition of new 
options. For example, a different profit function (eg. for a 
specific waferboard mill) could easily be incorporated into 
MAXPRESS simply by defining the new relationship in APL and 
then replacing the old function with this newly defined 
version. This inherent flexibility of APL was the primary 
reason for its selection as the computer language for the 


development of MAXPRESS. 


6.3 Model Components 
As mentioned previously, MAXPRESS is characterized by a 
modular composition. The main components, defined by their 


APL function or variable name, are described here. For a 
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complete, fully commented alphabetical listing of all APL 


functions in MAXPRESS, refer to APPENDIX II. 


6.3.1 Main Routine 

The principal component of MAXPRESS is a function named 
SEARCHMAX. This function performs the pattern search, prints 
intermediate and final output, and calls the other functions 
in MAXPRESS, as required. SEARCHMAX is called by the 


function START. 


6.3.2 Gap Search Routine 

GAPQ, a function called by SEARCHMAX, performs a_= gap 
search, if so desired. The user specifies the perturbation 
depth desired. This refers to the number of steps away from 
the last Lagrange solution that the user would like to 
consider (a maximum of three is allowed). For example, with 
a perturbation depth of one, the Gap Search Routine would 
examine all possible combinations of the decision variables 
(moisture content, resin content, panel density and press 
time) + one step away from the base set of values specified 
by the user. The optimal value (i.e. not dominated; see 
Method for Handling Gaps ) nearest to, and greater than, the 
minimum level of internal bond specified by the user is 
selected, and the appropriate output for this combination of 
variables is printed. The perturbation depth requested when 
the Gap Search Routine is invoked should depend on _ several 


Parameters. In general, the smaller the gap, or the closer 
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the desired level of internal bond is to the last solution 
level of internal bond, the smaller the perturbation depth 
required to produce the optimal solution. In addition, the 
perturbation depth requested should depend on the size of 
the computer used. Memory requirements, and the cost of 
running the Gap Beeeeh Routine, increase rapidly with 
increasing perturbation depth. Ultimately, the selection of 
perturbation depth is subjective in nature and depends upon 


past experience of the user. 


6.3.3 Objective Function 

The objective function (OBJFCN) used in  MAXPRESS 
computes the variable profit (dollars/8-hour shift) net of 
variable expenses. Fixed costs are not considered. These 
terms of reference are completely arbitrary. If desired, a 
user could easily re-define OBJFCN with different terms of 
reference. It should be noted, however, that a true fixed 
cost (i.e. one which must be borne regardless of the 
activity level of the waferboard mill) will not affect the 
solutions derived uSing MAXPRESS. To illustrate how a simple 
objective function can be constructed, the derivation of 
OBJFCN will be discussed. 


Definitions: 


tes TPS Total’ Profits 
2. VP ="Variable Profit, 
S.eTerazeTOtal=Costs- 
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5. VC = Variable Costs, 

6. TR = Total Revenue, 

7. VR = Variable Revenue, 

8. Q = Quantity produced; expressed as MSF/Shift, 

9. t = time for one press load to be completed (minutes), 


10. T = Time in one shift (minutes). 

The two basic components required in the objective 
function are expressions for VR and VC. A Suitable 
"activity' base must be chosen for the allocation of these 
costs and revenues. Fluctuations in the activity base should 
be closely correlated to fluctuations in the variable costs. 
and revenues. In OBJFCN, one MSF (thousand square feet; of 
whatever panel thickness one is interested in) was chosen 
for the activity base. Q must then be calculated (one could 
refine this calculation by allowing for down time; see 
CONCLUSION AND RECOMMENDATIONS) : 

Q (MSF/S) = T/t x 48 (Panels/Press Load) x 0.032 (MSF/Panel) 

Once Q has been established, VC and VR must _ be 
calculated. VR/MSF is provided by the user, and is simply 
the selling price/MSF at the factory gate for a particular 
thickness of waferboard. VC/MSF is somewhat more Beer cure 
to calculate. Following are the different components of VC, 
together with the decision variables (resin content = RC, 
moisture content = MC, press time = PT, and panel density = 
DEN; see Decision Variables - The Press Cycle) which 
directly affect these costs: 


ineuResin cost. =) :£( RC) «DEN)» 


93% @e@?. ~~ @@GGAI9Q¥,.. c0)mscra) 


tS oe gms ; 
see. best ——— 
: online 
Lope nena ea 


. wis 
; , 

py » 0 

20 ea 
uc / 99 patlems 

ry ~ 3 @& 

ise A '- 
a /@a)06e0, i i ie 


4 @6%4 a 14 @ dom 19D 


200624 esees ette < 
$085 .a)},4 


a 7 
_ wilte,. | 


48 


2. Wax cost = f£(DEN), 

Be Woodkocost ="f(0DEN), 

ae Fuel cost (forudrying waters) Y= Gf (MC; PDEN) 

5. Other variable costs = f(overall mill activity) 
ala DOUT? 
- waste (eg. panel trim, wafer fines, etc.), 
- variable overhead, | 
-~selectricrty, 
- other supplies, 
- miscellaneous. 

The decision variables included in the optimization 
model interact with each other to affect VR and VC, both 
directly and indirectly. Press time is the only decision 
variable which directly affects VR, but changes in the other 
variables will indirectly have an effect through their 
impact on panel quality. The effect on panel quality could 
change the optimal level of press time, hence, this could 
also change VR. Similarly, resin content, moisture content, 
and panel density have a direct effect on VC, but press time 
can indirectly have an effect through its impact on panel 
quality, and subsequent changes in the optimal levels of the 
other decision variables. 

VaeeOBJECN). (ontypathe a firstisioureyitems tiny the sist, 
above, are treated explicitly, as they can be easily 
allocated to MSF. The remainder are lumped into a category 
cCalled@picthergavariableyrcosts’. Splitting “this general 


category into smaller, more specific, categories could 
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result in a significant improvement of OBJFCN, if they could 
be allocated properly. VP can then be derived: 

Mei dollars/shift) =.O7xeiVR 7 VG} 
If FC is known, TP can be calculated: 


TP (adolwars/Shi ft) =) VPi—-FC 


6.3.4 Constraint Data Variable 

The constraint data variable (CONMAT) used in MAXPRESS 
was derived from various sources in the literature. Internal 
bond was selected for use as the only quality constraint in 
MAXPRESS. Although several measures of quality could be used 
in the Lagrange optimization of the press cycle, a different 
value for LAMBDA would have to be selected for each, 
complicating the operation of the model considerably. In 
addition, from an operational viewpoint, it appears that 
internal bond is usually the limiting factor in waferboard 
production. Attainment of adequate MOR and MOE is usually 
much easier than reaching adequate levels of internal bond 
(James 1981). 

The primary source for relationships between panel 
quality (internal bond) and production variables was Udvardy 
(1979) (for a graphical representation of this data, refer 
to APPENDIX III). The panels in this study had the following 
specifications: | 
- 7/16" aspen waferboard, 
= Landustry producedsefiakes], 


- press temperature = 400°F., 
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PenoLscune Content = "6, Of OFD? -Lurnish, 

mCesSTMeCONntCente= zy and 27 on rot wO.De furnish *phenol 
formaldehyde resin), 

- nominal panel density = 32, 40, and 48 pcf, 

- press time = 3, 5, 7 minutes. 

More data was sought, especially for different levels 
of production variables which were not varied here (eg. 
press temperature, moisture content, flake dimensions, 
etc.), but with little success. Difficulty was encountered 
due to a lack of published data and also to the wide variety 
of testing methods used by various researchers. A reasonable 
approximation for different levels of moisture content was 
obtained by averaging the results of several studies (Johns 
Shite ey ea O tr, ae gan and Schiewind 1971, Bryan and 
Schiewind 1971). The published data for the four resulting 
production variables (moisture content, resin content, panel 
density, and press time) were used as reference points 
through which relatively smooth curves were hand-fitted. 
Specific points (i.e. specific values for internal bond, 
measured in psi, for different combinations of the four 
Decision variables) mewere then enteredsintO aidatasvariable 
for use with MAXPRESS. The resulting data variable has_ the 
following dimensions: 5 (levels of panel moisture content; 
from 6% to 8% on an O.D. basis) x 5 (levels of panel resin 
content; ‘from 2% to 2.5% phenol formaldehyde resin) x 9 
(levels of panel density; from 32 to 48 pcf) x 9 (levels of 


press time; from 3 to 7 minutes) = 2025 elements. These data 


ie 


are only a rough approximation to what one might use ina 
real-world application of MAXPRESS, and are only intended to 
demonstrate the operation of the model. If better data were 
available (eg. from a waferboard production company), they 
Should be incorporated into the model. In addition, more 
production variables could be added if the data were 


available. 
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7. MODEL USER'S GUIDE 


7.1 Model Inputs 

Two types of input are required to operate MAXPRESS. 
The first type of input must be prepared beforehand, and 
loaded into an APL workspace with MAXPRESS. This type of 
input is referred to as Ceneral input mune seconditype of 
input is provided at the user when MAXPRESS is run. This 


type of input is referred to as 'User Provided Input'. 


7.1.1 General Input: 

There are two major pieces of general input required to 
run MAXPRESS; 1) an objective function and; 2) a constraint 
relationship (in the form of discrete data). The objective 
function provides the payoff relationship for the particular 
mill in question, while the constraint data provides the 
quality relationship between decision variables (eg. press 
time, resin content etc.) and some measure of quality such 
as internal bond (see Model Components for a description of 
the standard objective function and constraint variable 
developed for MAXPRESS). The objective function is defined 
as an APL function (called 'OBJFCN'; see APPENDIX II), while 


the constraint data is stored in an APL data variable. 
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7.1.2 User Provided Input: 


MAXPRESS prompts the user for the following cost and revenue 


inputs: 


1. Panel thickness (inches; currently, the user can only 
use §.4375" (7/16™)" due to a lack of data -for panel 
quality relationships to production variables), 

2. Resin cost (dollars/pound), 

3. Wax cost (dollars/pound), 

4. Wood cost (dollars/O.D. pound of wafers), 

5. Fuel cost (dollars/MCF natural gas), 

6. Other variable costs (dollars/MSF waferboard), 

(parse iangptornice sat) the.) Lactory, -qaatew (dollars/MSE 


waferboard), 


The minimum possible value, maximum possible value, and step 


size 


(for the Hooke-Jeeves Direct Search Algorithm - step 


size must equal the interval size of the constraint matrix), 


for each of the following variables: 


The 


@s Moisture content. of ‘the’ panel» (% ‘of O.D. panel 
weight), | 

9. Resin content of the panel (% of O.D. panel weight), 
10. Panel density (pcf), 

17s. Press time (minutes), 


Starting point for the Hooke-Jeeves Direct Search 


Algorithm: 


(ope Moisturercontenc (4) 
fSeResin econmtenta Gn), 


14. Panel density (pcf), 
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15. Press time (minutes), 

16. Minimum quality (eg. internal bond) desired, 
Parameters required for Everett's method of Lagrange 
Multipliers: 

17. The value for LAMBDA, 

Miscellaneous: 
18. YES/NO answers to questions posed by the model. eg: 
"Would you like to continue this analysis?' 


"Would you like to change the value of LAMBDA?'. 


7.2 Model Outputs 

Output from MAXPRESS includes production resource levels: 
1. Moisture Content (% of O.D. panel weight), 
2. Resin Content (% of O.D. panel weight), 
3. Panel Density (pcf), 
4; Press Time (minutes), 

Parameters relating to the Lagrange optimization: 
5. Value of the Lagrange Function, 
6. Value of the Lagrange Multiplier (A; printed at start 
and end of analysis only), 
We Value of the Objective (Payoff) Function 
(Dollars/8-hour shift; net of variable costs but does 
not include fixed costs), 
8. Value of the Constraint (Quality) Function (Internal 
Bond in psi for base case, but any other suitable 


measure of quality such as MOE or MOR could be easily 
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substituted). 
Output iSsapranted. at ~the.. start ‘Re the analysis, at 
intermediate solutions (base points) in the pattern search 
routine, at the optimum solution for the particular value of 


X used, and when the Gap Search Routine is invoked. 


7.3 The Role of LAMBDA 

LAMBDA (4), the Lagrange multiplier, represents the 
marginal value per unit of resource (panel quality, measured 
by internal bond; see Economic Significance of the Lagrange 
Multiplier in Chapter 5). Therefore, raising the value of 
LAMBDA will generally raise the level of internal bond in 
the optimum solution provided by MAXPRESS because the user 
has implied that panel quality has more value (hence, more 
quality and a higher level of internal bond at the optimum). 
The reverse is true when LAMBDA is lowered. An exception to 
this occurs when the range of values for LAMBDA crosses a 


gap region. 


7.4 Using the Model - A Sample Run 

In order to better demonstrate the use of MAXPRESS, the 
following narrative of a sample run of MAXPRESS is provided 
(the actual run is found in APPENDIX IV). 

After entering APL mode, and loading the appropriate 
program and data workspaces, the user enters the command 


"START', to initiate execution of MAXPRESS. The program 
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prompts the user for required inputs, and then prints 
initial, intermediate, and final output results. At this 
point, the user may stop the analysis, invoke the Gap Search 
Routine, or change user specified parameters and re-run the 
program. 

In the example provided (APPENDIX IV), the user enters 
the required input (pages 103 and 104), and receives output 
for this mix of inputs and value of LAMBDA (24) (page 105). 
The resulting value for internal bond is 32.909 psi, with 
variable profit/shift of $10150. As an example, consider the 
desired level of internal bond to be 42 psi (the CSA 
standard is 40.6 psi; National Standard of Canada 1978). To 
achieve this level, the model must be re-run, with a higher 
value for LAMBDA. In the example, LAMBDA is raised from 150 
to 160. The solution (page 106) shows internal bond of 
48.282 psi, and variable profit/shift of $7810. This level 
of internal bond is somewhat high, so the model is re-run 
with LAMBDA of 155 (page 107) because this implies a lower 
value for panel quality. The resulting solution is the same 
as for the previous value of Lambda, so the model is run 
again with LAMBDA of 152.5 (page 108). The solution for this 
level of LAMBDA has internal bond of 41.068 psi, and 
variable profit/shift of $8910. This is very close to the 
desired level of LAMBDA, but if 42 psi was considered the 
absolute minimum level of internal bond allowed, the model 
would have to be re-run. The value of LAMBDA could be 


adjusted upwards Slightly, and the model re-run as before, 
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but this is a rather time consuming effort considering how 
close the last level of internal bond is to the desired 
level. Instead of relying on trial and error to adjust the 
level of internal bond by such a small amount, the Gap 
Search Routine was invoked using a perturbation depth of one 
(page 109). The resulting solution has internal bond equal 
Wome 07 6. .0S1 wancevarhable prota tyshit weotes coves) This 1s 
the optimal solution for this selection of input parameters, 
and is as close to the desired level of internal bond as can 
be obtained with the discrete constraint data used (see 
Method for Handling Gaps in Chapter 5). This search (and the 
searches in the remainder of this analysis) was started from 
several different base points, once the level of internal 
bond was close to the desired level, to help ensure that a 
true global optimum was found. Because the global optimum 
had been found, in all cases, output from these searches is 
not included in APPENDIX IV. 

This analysis was continued (page 110), with the cost 
of resin being raised from $0.70/pound (first case) to 
$1.00/pound. The value of LAMBDA remained the same (152.5), 
and the model was re-run (page 111). In this case, internal 
bond was 45.892, and variable profit/shift was $6487. Since 
this level of internal bond was close to the desired level 
of 42 psi, it was again decided to use the Gap Search 
Routine, this time using a perturbation depth of three, 
instead of one as in the last case, because of the greater. 
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psi) and the desired level of internal bond (42 psi) (page 
i112). The »~resulting™ solution ‘provides "internal ** bond “of 
ee e0Se «psi, and” variable »profit/shitt ™ of #s6762. "lt is 
interesting to note that a jump in resin cost of $0.30/pound 
(from $0.70 to $1.00) resulted in a decrease of variable 

Brorit/Snift of VS 1814 (from $8576 to $6762), and changed the 
cphenneal levels of moisture content (from 7% to 7.5%), resin 
eoncent  (from?2.5%ato 2.125%), "and press eanemceeom AO PATL Ty 
to 4.5 min.). This illustrates the sensitivity of optimal 
operating policy (and profits) to changes in the costs of 
resin, one of the four decision variables. 

This run wasS continued with changes to the input 
Parameters (page 113), in order to illustrate a gap. The 
first solution, with LAMBDA at 90 (page 114) had internal 
bond of 755-598 Spsi “and variable @ *profit/shite of $1655: 
LAMBDA was dropped to 87, and the model re-run. This 
produced a solution (page 115) with internal bond of 1.800 
psi and variable profit/shift of $7490. Another run, with 
LAMBDA at 88, produced a solution (page 116) that was 
PTAOeALICcal. “‘Co=rthe. sfirsteyesolution. Pwithe -GAMBDAY at 90. 
‘Subsequent runs (not shown in APPENDIX IV), with the value 
of LAMBDA varied between 87 and 88, failed to produce any 
new base points. This indicated that a gap existed in the 
area of interest, and that further manipulation of LAMBDA 
would not likely, produce “an’*acceptable™ solution. =)In 
practice, the author found that reducing the difference 


between consecutive levels of LAMBDA to an interval of less 
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than one unit (eg. between 87 and 88 in the case above) was 
not justified because of the effort this requires, and 
because, in all cases tested, no new base points were 
discovered using this strategy. Because the gap was large, 
another run was made (page 117) to ensure that no base 
points existed which would produce a value for internal bond 
closer to the desired value. Minimum internal bond, a 
Parameter which should be set at 0O for standard runs of 
MAXPRESS, was set at 42. When the pattern Search reaches the 
minimum level for internal bond, it performs exploratory 
moves only. This action could, in some instances, reveal a 
new base point. It is conceivable, with some irregular 
functions, that a true base point (i.e. the true maximum of 
ene Lagrange function for a particular value of LAMBDA) 
could be skipped over by a pattern move, particularly if the 
pattern move is large. Pattern moves get larger, or 
‘accelerate’, once a pattern is established (see Pattern 
Moves), but exploratory moves are always the same size 
(equal to the step size of the discrete constraint data) in 
MAXPRESS. Thus, it is impossible for an exploratory move to 
skip over a true basepoint. A true base point would have a 
higher value for the Lagrange function than the previous 
solution obtained through a normal application of MAXPRESS. 
Finding a new base point would narrow the gap, thereby 
helping to find the best solution for the quality level 
desired, even if the Gap Search Routine still has to _ be 


used. In the example, no additional base points were 
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discovered. This strategy would have to be used only where 
gaps are very large (a subjective assessment, which can only 
be made on the basis of experience with MAXPRESS). In fact, 
the author discovered only one case, in numerous runs of 
MAXPRESS, where a new base point was revealed using this 
strategy. In the example, the gap was small enough to 
evaluate directly with the Gap Search Routine, and this was 
Subsequently done (page 118). The final solution has 
imternat bondmofe+2.1t8tandmvariablesprotirt Zshut uno af2855% 


At this point, program execution was terminated. 
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8. MODEL TESTING 

Two components of MAXPRESS were tested; the main 
routine, and the Gap Search Routine. Testing of the main 
routine involved the arbitrary selection of a differentiable 
function with three independent variables for use as the 
objective function. Another differentiable ee iont also 
with three independent variables, was chosen as_ the 
eonstraint «Eunctlon™ and was used toil! aes0 x9 50)x«"30 
(= 45000 elements) data matrix (one dimension for each 
independent variable), using an arbitrarily selected step 
Size for each independent variable. MAXPRESS was then used 
to optimize the objective function, subject to the 
constraint data matrix. These results were compared to the 
results obtained by solving for the optimum using 
differential calculus. MAXPRESS provided results comparable 
to those obtained directly. Naturally, there was a= slight 
difference due to the discrete nature of the constraints, 
and the subsequent discrete nature of solutions. 

Several gapS were encountered using the objective 
function and constraint data variable developed for this 
project. Testing was conducted on one of these gaps by 
constructing a complete list of possible solutions, 
selecting the optimal solution, and comparing this to the 
result provided by the Gap Search Routine. For the _ gaps 
tested, the Gap Search Routine always provided the optimal 
solution with a perturbation depth of three and, in most 


cases, a perturbation depth of only one or two was required. 
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MAXPRESS appears to perform well, solving the trial 
problems in all cases tested, but it is not absolutely fail 
proof. While the model performed well in all of the tests 
and trial runs conducted, it is possible that non-optimal 
solutions could be generated, particularly with the use of 
unusual constraint and objective relationships. However, if 
reasonable precautions are taken (Such as starting the 
pattern search from various base points, once it is thought 
that an optimal solution has been found), the generation of 
true optimal oleh one should almost always occur. 

A final observation is that MAXPRESS is easy to use. 
Once the initial objective function and constraint data 
variable are set up (and this would have to be done just 
once for each mill configuration), the user has only to 
respond to queries from the computer, and vary the level of 
LAMBDA to achieve the optimal solution for the desired level 


of internal bond. 


i a a a? Lhd 
) 


felis? ets ealeleas 


ites yisdviceds. rus 
eiae2. Geis ae t io 


oe "si be a. Aizawa sy 


G 4 seotves 2h 
p59 Con el iwine 
‘ sete en 
bake | ‘-<élagesniee 


: ser oserd 


Pr er + a 


9. CONCLUSION AND RECOMMENDATIONS 
Based on the testing conducted, and the numerous trial 
runs undertaken with various inputs, it appears’ that 
MAXPRESS performs very well in the constrained optimization 
of the waferboard press cycle. 
In addition, MAXPRESS meets four specific design 
considerations established ae the inutiationes Of ethas 


project. Namely: 


MAXPRESS accepts non-differentiable, non-linear, and 


discontinuous objective functions; 


MAXPRESS accepts constraint relationships in the form of 


discrete data; 


MAXPRESS is an optimization, not simulation, model; 


MAXPRESS is easy to use, even for personnel with no formal 

training in computer programming or operations research. 

As it exists, MAXPRESS has several potential uses, but 
these will depend on a refinement of both the objective 
function and the constraint data variable. The usefullness 
of the procedures developed for the optimization of the 
press cycle is far more important than the actual profit 
function or panel quality data assembled for this project. 
One major use could be the identification of optimal 
operating —policyy sini stermseyok resinetcontent.{ moisture 
content, press time, and panel density) ain response to 
changeseeinescostiiton diiierentatactorSyomionrodicti onwabon 
example, one might pose the question; ‘If resin prices 


increase 15%, should operating policy be changed (perhaps 
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through lower resin content, and higher press time)?' 

Another significant use could be the testing of different 

resins based on their price, and on what is known about 

their properties. Other sensitivity analyses could be 
conducted to test the effects of changes in other factors of 
production. 

Several recommendations can _ be made LOr future 
research: 

1. A complete waferboard mill model should be constructed, 
using MAXPRESS as a central portion, so that the entire 
production process can be examined. Mill design and 
re-design questions could then be addressed, and the 
implications of operating policy decisions resulting 
from MAXPRESS could be assessed from the entire mill 
point of view. For instance, changes in press time will 
affect the entire mill throughput. A complete mill model 
would allow one to gauge the impact of these changes, 
and to predict whether the rest of the mill could keep 
up with the press. 

The complete mill model should probably be a 
queuing model based on the work ae Carino ‘and ‘Bowyer 
(1979, 1981). Their model provides optimization, using 
the Hooke-Jeeves Direct Search Algorithm, as well as 
Simulation of the queuing system. The work of Rosenshine 
and Chandra (1975a, 1975b) could be useful for queues 
which have batch arrivals (such as the press). 
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(payoff) function and the constraint data variable, as 
they are currently only preliminary versions designed to 
illustrate the operation of MAXPRESS. If other 
production? variables» jare deemedzeimportant to panel 
quality (eg. press closure rate and pressure, wafer 
dimensions etc.), appropriate data should be assembled 
and included as new dimensions to the constraint data 
variable. In addition, data should be assembled for 
different panel thicknesses. 
Refinements EO the objective function could be 
implemented to allow the comparison of various fuels (or 
wood waste asa fuel). In addition, it might be useful 
to break down the ‘other variable costs’ category into 
several components, and allocate them 2 eee to 
production. Another significant improvement would be to 
allow for variable down time, and variable waste 
factors. The waste factors could represent losses at 
different stages of the waferboard production process 
(eg. loss of panel at the trim saws; loss of wood at the 
debarker; etc.). 

phewerchangessssuggestedsyheres wouldem probably — be best 
achieved by working in cooperation with personnel of a 
waferboard production company. | 

New panel quality constraints should be added to 
MAXPRESS, or the current constraint (internal bond) 
should be replaced, if other constraints are deemed 


important. Other useful panel quality constraints could 
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include MOR and MOE. 

4, MAXPRESS should be translated into FORTRAN and, possibly, 
other popular languages as well. While APL is extremely 
useful for developing an optimization model, it is not 
the best language to promote to waferboard production 
companies. APL is not available on many small computers, 
and where it is available, the workspace size (APL 
equivalent of memory), would likely be inadequate to run 
MAXPRESS. A FORTRAN version of MAXPRESS might also _ be 
cheaper to run than the APL version. 

5. A useful option to include in MAXPRESS would be the 
possibility for a complete breakdown of production costs 
into categories such as fuel, resin, labour etc. This 
option would require an improved objective function. 

In conclusion, it 1S evident that MAXPRESS currently 
has several potential uses, and that there are several areas 
where the model could be significantly improved. Throughout 
this study, emphasis has been placed on producing a 
practical, usable tool for the waferboard industry. MAXPRESS 
appears to be such a tool, but its ultimate test will be 


acceptance and use by the waferboard industry. 
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12. APPENDIX II - LISTING OF FUNCTIONS USED IN MAXPRESS 
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IAI E Daet les 2S thse Obi. 2131/ 62 


- MAXPRESS 
Pe AUG 44,1982 
8:54:58 
xxxx GROUPS xxxxxx 
DIRECTORY FOR MAXPRESS 
RHEE OE YU ON COT LION (8 CRESS 
BASEQ GAP LTR PARAM 
BY GAPQ MATRIX1 ROW 
COL TF MATRIX2 SEARCHMAX 
CONFCN IFNOTO MATRIX3 SELECT 
COSTQ IFYESTO MESSAGE START 
EXPLAIN INITIAL MINQ STATUSQ 
EXPLOREMAX INTERPOLATE OBJFCN STEPSET 
FUNCTION LAMBDAQ ON1 VARIABLES 


ease VY AOR POAC BILE Sere & ex 


( NAMES WITH UNDERLINED CHARACTERS ARE NOT LISTED.) 
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B44? Sb 406 © 


+. tae 
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xxxx  BASEQ BASEQ 
0| VY BASEQ 
Ue 
2! +0 IFNOTO ‘WOULD YOU LIKE TO CHANGE THE STARTING POINT FOR THE SEARCH? 
3| FIRSTBASE<INITIAL, 1+FIRSTBASE 
yl om 
S| a 


6| a ASKS QUESTION FOR CHANGING STARTING POINT OF SEARCH 


wi 
kK KK BY BY 
5 BY CCCI ON xpa)iiipa) tai, (Gi Oxpiainpavt a 
2 :(0=0\0/,a+COL a)#(0=0\0/,w<COLw) 
3 | 3(¥a) BY (¥a) 
4y| oa 
5| a CREATE A MATRIX FROM a AND w 
6| a BY FIRST CONVERTING THEM TO MATRICES AND THEN 
7| a ADJUSTING THEIR ROW SIZES TO MATCH 
8| a AND PLACING a TO LEFT OF ao. 
9| mn IF THEY DIFFER IN TYPE THEY ARE CONVERTED 
101 a TO CHARACTER FORM. 
xxxx COL COL 


(KVES C26 (IA) Sak A A) 

A 

a RESTRUCTURE w AS MATRIX WITH AT LEAST ONE COLUMN. 
A (ONLY FIRST TWO COORDINATES OF STRUCTURE ARE 

a RETAINED.) 
A 
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aT 
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- a - 7 
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xxxx  CONFCN CONFCN 


VY OUT<CONFCN B;COUNT ;MAX;:COUNT:I3MAXI3MATR3MAX;3COUNTER:3V3W 
COUNTER+<1 

B<ROW B 

MAXI<1itoB 

MAT R<(1+pB)p0 

BEGIN: 

>END IF COUNTER>MAXI 

V<B[COUNTER; } 

I<(p,V)p0 

MAX<p,V 

COUNT<1 

dal |} She ;079 

12| +SET IF COUNT>MAX 

13| ILCOUNT)<(V[COUNT ]=CONSTEP[COUNT;1)/.1pCONSTEP[COUNT; J 
14| COUNT<COUNT+1 

15| +START 

NGaleeSiees 

17| MATR(LCOUNTER]<CONMAT(I(11];I2(02]3;7(31];7(4]] 

18| COUNTER+COUNTER+1 


OWOWDAYMMNGCWNF OC 


19| +BEGIN 
20| END: 
21| OUT<MATR 
22 A 
23| a CONFCN SELECTS THE APPROPRIATE VALUE FOR INTERNAL BOND 
24| a FROM THE CONSTRAINT MATRIX (CONMAT) BY LOCATING EACH 
25| a OF THE FOUR DECISION VARIABLES (RC, MC, DEN, PRESS) 
26! a ALONG THE CONSTRAINT STEP SIZE VECTORS FOUND IN 
27| a CONSTEP. 
DiStieA 
xxxx COSTQ COSTQ 
0| V COSTQ 
aah 3a 
2| 270 IFNOTO ‘WOULD YOU LIKE TO CHANGE THE VARIABLE COST FIGURES?' 
3| COST<PARAM 
4 | ey. 
5l a 
6| nan ASKS QUESTION FOR CHANGING THE VARIABLE COST FIGURES. 
valle . 
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"DESIGNED FOR THE OPTIMIZATION OF THE PRESS CYCLE' 
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‘TWO TYPES OF INPUT ARE REQUIRED TO RUN MAXPRESS:' 

‘1. GENERAL INPUT 

‘2. USER PROVIDED INPUT." 
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‘QUALITY RELATIONSHIP BETWEEN DECISION VARIABLES (EG. PRESS' 
‘TIME, RESIN CONTENT, MOISTURE CONTENT, PANEL DENSITY) AND! 
‘SOME MEASURE OF QUALITY SUCH AS INTERNAL BOND. THE OBJECTIVE' 
‘FUNCTION IS DEFINED AS AN APL FUNCTION CALLED OBJFCN, WHILE' 
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en 
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USE .4375 INCHES (7/16) DUE TO A LACK OF DATA FOR PANEL’ 
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' 9. RESIN CONTENT OF THE PANEL (PERCENT OF 0.D. PANEL WEIGHT) ,' 
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‘THE STARTING POINT FOR THE HOOKE-JEEVES DIRECT SEARCH ALGORITHM: ' 
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MOISTURE CONTENT (PERCENT),' 

RESIN CONTENT (PERCENT),' 

PANEL DENSITY (PCF),' 

PRESS TIME (MINUTES) ,' 

MINIMUM QUALITY (EG. INTERNAL BOND) DESIRED; SHOULD BE SET' 
TO 0 IN MOST CASES; ONLY EXCEPTION IS WHERE YOU WANT TO' 
EXPLORE THE AREA OF INTEREST USING EXPLORATORY MOVES TO' 
SEE IF THERE ARE NEW BASE POINTS; FOR INSTANCE, IN GAP ! 
REGIONS,' 


PARAMETERS REQUIRED FOR EVERETTS METHOD OF LAGRANGE MULTIPIERS: ' 


‘17. THE VALUE FOR LAMBDA, THE LAGRANGE MULTIPLIER,' 

tha) 

"MISCELLANEOUS: ' 

"18. YES/NO ANSWERS TO QUESTIONS POSED BY THE MODEL.' 

Ua 
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EXPLOREMAX EXPLOREMAX 


V EXPLOREMAX;COUNT ;M 

COUNT<1 

M+p (~1+7EMPBASE ) fe 
START: 

>0 IF COUNT>M 

TEMPBASE( COUNT J<TEMPBASE( COUNT J)+STEPSIZE(COUNT ] 
>OUT IF( FUNCTION TEMPBASE)>TEST 

TEMPBASE( COUNT ]<TEMPBASE( COUNT ]-2xSTEPSIZE(COUNT ] 
>OUT IF(FUNCTION TEMPBASE)>TEST 

TEMPBASE(COUNT }«TEMPBASE([COUNT ]+STEPSIZE(COUNT } 
>COUNTER 

OUT: 

TEST«FUNCTION TEMPBASE 

COUNTER: 

COUNT<COUNT+1 

>START 

a 

EXPLOREMAX PERFORMS EXPLORATORY MOVES FOR THE 
HOOKE-JEEVES DIRECT SEARCH ALGORITHM, WHERE THE 
OBJECTIVE IS FUNCTION MAXIMIZATION. 


DoDD 


FUNCTION FUNCTION 


V OUT<FUNCTION V3;L;COUNT 

L<~14V 

Vi<merley iV 

>RESET IF(v/V<MIN)v(v/V>MAX) 

a >INTER IF MODE=1 

OUT<(COST OBJFCN V)-Lx(CONFCN V) 

>0 

INTER: 

OUT+(COST OBJFCN V)-Lx(INTERPOLATE V) 

70 

RESET: 

OUT<TEST 

A 

FUNCTION CHECKS THE DECISION VARIABLES TO SEE IF 
THEY ARE WITHIN THE ACCEPTABLE RANGE, AND THEN 
CALCULATES THE LAGRANGE SOLUTION FOR THIS COMBINATION 
OF DECISION VARIABLES. 

OBJFCN IS CALLED TO CALCULATE THE PROFIT/SHIFT. 
CONFCN IS CALLED TO DETERMINE INTERNAL BOND. 
INTERPOLATE IS NOT CURRENTLY USED, BUT COULD 

BE VALUABLE IS LINEAR INTERPOLATION BETWEEN 
STEPS OF THE DECISION VARIABLES WAS DEEMED IMPORTANT. 
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VY OUT<B GAP Q;BASE1;0B;CO;L;DIFF;BASE 
BASE<B 

BASE1+(COST OBJFCN BASE)-LAMBDAXCONFCN BASE 
OB<+COST OBJFCN Q 

CO<CONFCN Q 

L+OB-LAMBDAxCO 

DIFF+L-BASE1 

QUT -—DIEFSBY SES BY OB SBYOCOSBY 70 
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PERFORMS THE SAME TASKS AS FUNCTION, 
BUT FOR THE GAP SEARCH ROUTINE. 
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>END IFNOTO ‘WOULD YOU LIKE TO PERFORM A GAP SEARCH?' 
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BASE<INITIAL 
MINQ 0 
LAMBDAQ 0 


Wak 


‘ENTER THE VALUE FOR PERTURBATION DEPTH (MAXIMUM IS 3)' 


NO 

>ONE IF N=1 

>TWO IF N=2 

*+THREE IF N=3 

ONE: 

A<BASE SELECT BASE GAP MATRIX1 BASE 
>CONTINUE 

TWO: 

A<BASE SELECT BASE GAP MATRIX2 BASE 
>CONTINUE 

THREE: 

A+BASE SELECT BASE GAP MATRIX3 BASE 
+CONTINUE 


CONTINUE: 
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BASEPOINT<A[5 6 7 8],LAMBDA 


OUT<1 


a THIS IS THE GAP SEARCH ROUTINE. 
aA APPROPRIATE DATA IS OBTAINED FROM THE 
a USER. THEN, FUNCTIONS ARE CALLED DEPENDING 


a ON THE PERTURBATION DEPTH REQUESTED. 
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IBS IBY 


IF:w/a 
A 


A RETURNS a IF w IS TRUE OTHERWISE RETURNS EMPTY VECTOR 
A 


IFNOTO IFNOTO 
V R<«LABEL IFNOTO QUERY;ANS;T 

MICOUE RY 0) me 

ANS<4+(ANS#' ')/ANS<1) 

>(R<(4/(ANS='N '))va/(ANS='NO '))/FIN 

>((A/(ANS='Y "))va/(ANS='YES '))/FIN 


>(A/ANS='STOP')/STP 

‘IMPROPER RESPONSE, REPLY Y OR YES FOR YES, N OR NO FOR NO' 
‘TRY AGAIN' 

<7 

Suh. OF! 

"EXECUTION INTERRUPTED, REMEMBER TO CLEAR THE STATE INDICATOR' 
>0 

FIN: R<R/LABEL 

A 


a ASKS A QUESTION (QUERY), AND EXECUTES 

a A RESPONSE IF THE ANSWER (PROVIDED BY 

a THE USER) IS NO. 

A 

IFYESTO IFYESTO 
VY R«LABEL IFYESTO QUERY:ANS;T 

[KQUERY,' ' 

ANS<4+(ANSt#' ')/ANS<{) 

>(R<(A/(ANS=°Y '))va/(ANS='YES '))/FIN 

>((A/(ANS='N '))va/(ANS='NO '))/FIN 


>(A/ANS='STOP')/STP 

‘IMPROPER RESPONSE, REPLY Y OR YES FOR YES, N OR NO FOR NO‘ 
'TRY AGAIN' 

>2 

ENGR TR Sy oe 

"EXECUTION INTERRUPTED, REMEMBER TO CLEAR THE STATE INDICATOR' 
70 

FIN: R<R/LABEL 

A 

ASKS A QUESTION (QUERY), AND EXECUTES 

A RESPONSE IF THE ANSWER (PROVIDED BY 

THE USER) IS YES. 
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O| Y OUT<INITIAL;A 

1| ¢ 

2| ‘ENTER, IN ORDER, VALUES FOR EACH PROCESS VARIABLE:" 
3] V4] 

4| "MOISTURE CONTENT (PERCENT )' 

S| vr 

6| AO 

Fal t 

8| "RESIN CONTENT (PERCENT )' 

3 | ee 

10| A+A,O 

kali eg 

LAD) MAME JOBING (IXGI9)) V 

ihe} t) .% 

14| A<«A,O 

TS s les 

16| ‘PRESS TIME (MINUTES)' 

ih 7 AI ? 

18| A<A,0 

19| Cat 

20| OUT<A 

21! A 

22! a I/O FUNCTION FOR INPUTTING VALUES FOR 
23| a THE PROCESS (DECISION) VARIABLES. 
24| aA 
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INTERPOLATE INTERPOLATE 


V OUT<INTERPOLATE V;LOWER;UPPER;L;L1;U;U13;MAX;COUNT ; DIFF ;MARK 
L<U<L1<UV1<(pV)p0 

MARK<+0 

MAX<poV 

COUNT<1 

START: 

>SET IF COUNT>MAX 

>RESTART IF(+/V(COUNT]=CONSTEP[COUNT;])>0 

U1 (COUNT }<1+ (VL COUNT ]}<CONSTEP[COUNT;])/CONSTEP(COUNT;] 
ULCOUNT }+(U1 [COUNT 1=CONSTEP[COUNT;1)/.1pCONSTEP(COUNT ;] 
LUCOUNT }<U[COUNT]-1 

Lil(COUNT ]<CONSTEP(COUNT;L(COUNT))} 

MARK+COUNT 

COUNT<COUNT+1 

>START 

RESTART: 

UC(COUNT J<L[I COUNT 1<( VL COUNT 1=CONSTEP([ COUNT; ])/1pCONSTEP[COUNT; } 
U1 (COUNT ]+L1(COUNT }<0 

COUNT<COUNT+1 

>START 

SYIIER 

+SET1 IF MARK=0 

LOWER<CONMAT(L(1};L(02];L(03)3;204)] 
UPPER<CONMAT(U[1];U[2];U[3];U[4]) 

DIFF+UPPER-LOWER 
OUT+LOWER+DIFFx(V(MARK]-L1[MARK])*+U1[MARK]-L1i(MARK] 

>0 

SEaelae 

QUT<CONMAT WE elms (G2 Viste S lee eel a) 

A 

a CURRENTLY UNUSED. PROVIDES INTERPOLATION 

a CAPABILITIES FOR THE HOOKE-JEEVES ALGORITHM. 

A 


LAMBDAQ LAMBDAQ 


V LAMBDAQ B;L 

+SET IF B=0 

ft 

+SET IFYESTO ‘WOULD YOU LIKE TO CHANGE THE VALUE OF LAMBDA?" 
>0 


‘ENTER THE VALUE FOR LAMBDA’ 

rf 

LAMBDA<) 

L+0-LAMBDA 

FIRSTBASE<«(~1+FIRSTBASE) ,L 

La | 

a 

I/O FUNCTION FOR CHANGING THE VALUE 
OF LAMBDA, THE LAGRANGE MULTIPLIER. 


RY D 
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LTR: a+.Xw 


a LINEAR TRANSFORMATION, OR PRODUCT SUM OF a WITH w 


MATRIX1 


V OUT<MATRIX1 NsA3B3C3sD3sM3S3I 

are} 

S<STEPSIZE 
A<81p(27pN[4)+S04]),(27pN04)]),(27pNC4I)-SC4)) 
B<+81p(9pNL31+S5(3]),(9pN[3]),(9pNL3]-S[3]) 
C<+81p(3pN[2]+5([2]),(3pNL2]),(3pNL2]-S[2]) 
D<8ip(NC1)+S01)) N01), (N01I-S11]) 

McD eB Ya Gupb Yeeb arb YaaA 

LOOP: 

+END IF I> "1tpM 

W<( (MC 3IT]SMAXCLIJ)A(MC3IJ2MIN[I]))/i1itpM 
M<M(W; ] 

WEST 

+LOOP 

END: 

OUT<M 

A 

SETS UP MATRIX WITH ALL POSSIBLE COMBINATIONS 
OF DECISION VARIABLES FOR PERTURBATION 
DEPTH OF ONE. 


RpDDvDD 
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MATRIX2 MATRIX2 
V OUT<MATRIX2 N3A3SB3C3D3M3W3I 

eal oie 

S«STEPSIZE 
A+6259(125pNC4J+2xS[4)),(125pNL4]+S504]),(125pN04J),(125pN04IJ-ST4)), (12 
SpNC4J-2xS[4]) 
B<625p(25pN(3J+2xS[3]),(25pN(3]+S5(3]),(25pN03]),(25pNL3]-SL3]),(25pNL3 
LEA One EEC CA IO Gee ee 
A cieareeceane mieeroean iat. Gicaeael waliecmaa os 

M_DS BY CeBYerB eb YaA 

LOOP: 

>END IF I> "1tpM 

W<( (MU 3 IT]SMAX(IJ)A(M[>3;IJ2MIN[IJ))/11+pM 

M<M([W; ] 

Tie Tid 


>LOOP 


a SETS UP MATRIX WITH ALL POSSIBLE COMBINATIONS 
a OF DECISION VARIABLES FOR PERTURBATION 


a DEPTH OF TWO. 


4,1982 / 8:54:58 / «xxx MAXPRESS xxx PAGE 12 


a. 


ie 
oe at -= : 

. ve : 
| Yee See 


ae eo eu . | 


a » _ | cn 
_ - ue : 7 LP : 
= by —— lias a. 


yi S 
+ yeethihs ores 7 ee 


ai 2 eythl, jextoue 
207 ¢ ere ue Ober 1 te ate 
eo a 
= jp agreenges 1m i ; ie 
Olax es ae ane 
ieitene 1k 
erties ¢ 
syiae le, 
NS 48 
aL. 
rina dee 8 
io | 1H 0 Bee 
Ty ee 
4 tae. at 
£6 eM yee hd ele api. 


87 


xxxx MATRIX3 MATRIX3 
Ol V OUT<MATRIX3 NsA3B3C3D3M3W3I 
SU] SESH OHA RAe 
A\| tera 
31 A#2401p(343pN[4]+3xS(4]), (343pN041+2xS04]),(343pNT4]+S04)),(343pN04]), 
(343pNL4I)-S(4]),(343pNL4)]-2xS04]),(343pN[4]-3xS(4]) 
4] BH2401p(49pN(31+3xS(3]), (49pN(3]+2xS(3]),(49pNL3I]+S(13]),(49pNT3]),(49~ 
N(3]-S[3]),(49pNL3]-2xS[3]), (49pN(3]-3xS([3]) 
Bi] GEA CEA NORE I) 5 (CED Neel ID) a Hone ESE s CHM) 5 CLA ale 
SIEZAD) 6 (POME2YV SRS SIC 2a) CAL MOSEL 1) 
6| D<2401p(N(11]+3xS(1]),(N01]+2xS501]),(N011]+5[1]),(N01]),(N01]-S(1]),(N[1 
12x S) (ts ay ON) el 13 2c Siliel ly) 
7| M<D BY C BY B BY A 
8| LOOP: 
g9| +END IF I>~1itpM 
10! W+( (MU 3IT]<SMAXCLIJ)*(ML3IJ2MIN[IJ]))/.i1+ pM 


11| M+MCLW3] 


ste2i) Becszecsal 
13| +LOOP 
14| END 
15| OUT<M 
16| aA 


17| a SETS UP MATRIX WITH ALL POSSIBLE COMBINATIONS 
18| a OF DECISION VARIABLES FOR PERTURBATION 


19| an DEPTH OF THREE. 
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20) 
Dh} 
D|| 
23| 
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MESSAGE MESSAGE 


V MESSAGE 


tel 
Rt 
VK KK KKK KKK KK KK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK! 


Vk KK KK ROKK ROKK ko kok Kk kok KO Ok Kk kk kk ek kk KKK KKK KK KKK! 


'xx WAFERBOARD PRODUCTION OPTIMIZATION MODEL xx' 


"xx ALBERTA RESEARCH COUNCIL xx! 
tx FOREST PRODUCTS PROGRAM: FP-19 xx! 
tx x AUTHOR: TOM I. GRABOWSKI xx! 
tex LAST REVISION: JUNE 1, 1982 xx! 


VK KK KK KKK KK KKK KKK KKK KKK KKK KKK KKK KK KKK KK KKK KK KKK! 


VK KKK KKK KKK KKK KK KKK KKK KK KKK KKK KKKK KKK KK KKK KKK! 
' , 


’ J 


'xxx*x*WARNING: YOUR DATA FILE MUST BE COPIED INTO THIS WORKSPACE.' 


S8: 


‘ENTER ONE OF THE FOLLOWING KEYWORDS:' 


,¢ 


"EXPLAIN (GIVES A DESCRIPTION OF THE PROGRAM)' 


‘START (INITIATE RUNNING OF THE PROGRAM)' 


CIS 


A 

aA MESSAGE WHICH APPEARS ON THE SCREEN 

a (AS A LATENT EXPRESSION) WHEN THE USER 

a LOADS THIS WORKSPACE. 

A 

MINQ MINQ 
V MINQ B 


cP ae © aay <A) 


+0 IFNOTO ‘WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND?' 
SET: 


f 


"ENTER VALUE FOR MINIMUM INTERNAL BOND' 


tf 


MINIMUM+(] 


sf 


A 

a I/O FUNCTION TO CHANGE THE MINIMUM LEVEL OF 
a PANEL QUALITY (INTERNAL BOND). 
A 
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rFPOW ONAN FWNP 


BPR 


AUG 


OBJFCN OBJFCN 


VY OUT<A OBJFCN B;ODWEIGHT;MC;VARCOST;REV3EXP 
B<+ROW B 
ODWEIGHT+<1.25x(0.01x100-B(31]+BL3;2J]+1)xBL33]x32xAl1]412 
MC+(1.15x1892xA[5]#1001132)-(0.05x1892x(A[5]#1001132)xBL31]) 
VARCOST<(A([61]1x0.032)+ODWEIGHTx(A[2]xBL32]4+100)+(AL3]x0.01)+MC+AL4 ] 
REV<(4804B[34])x24x2x0.032xA(7] 
EXP<(480+B[34])x24x2xVARCOST 
OUT+REV-EXP 
a INPUT A: 1) PANEL THICKNESS (INCHES), 
2) RESIN COST (DOLLARS/POUND), 
3) WAX COST (DOLLARS/POUND), 
4) WOOD COST (DOLLARS/O.D. POUND CHIPS), 
5) FUEL COST (DOLLARS/MCF NATURAL GAS), 
6) OTHER VARIABLE COSTS (DOLLARS/MSF), 
7) SELLING PRICE (DOLLARS/MSF). 
INPUT B: 1) MOISTURE CONTENT OF PANEL (PERCENT OF 0.D. WEIGHT), 
RESIN CONTENT OF PANEL (PERCENT OF 0.D. WEIGHT), 
3) DENSITY OF PANEL (PCF), 
4) PRESS TIME (MINUTES). 


VARIABLE COSTS AND QUANTITIES OF PROCESS VARIABLES ARE INPUTS, 
VARIABLE REVENUE PER 8-HOUR SHIFT IS OUTPUT. 
ODWEIGHT ASSUMES 25 PERCENT OF WAFERS LOST TO TRIM, ETC. 


PrxrnvrrvrvrerDIBRePrex Rae re VRv,PMNyD 
N 
~~ 


ON1 ON1 


ON1 :(((0 1 xpw)[pa)ta), (D70)((0 1xpa)[ pw) tw 
:(0=0\0/,a<ROW a)2#(0=0\0/,w<ROWw) 
:(¥a) ONI (¥w) 


CREATE A MATRIX FROM a AND w 

BY FIRST CONVERTING THEM TO MATRICES AND THEN 
ADJUSTING THEIR COLUMN SIZES TO MATCH 

AND PLACING a ON TOP OF w. 

IF THEY DIFFER IN TYPE THEY ARE CONVERTED 

TO CHARACTER FORM. 


parvo vpvDRPrPDMDDDD 
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16 | 
9) || 
18 | 
19] 
20| 
2 
D2| 
23 | 
24 | 
25| 
26 | 
27 | 
28 | 
29 | 
30] 
shah || 
<1) | 
33] 
34 | 
35 | 
36 | 


KkKK*K 


mWMFwWNe 


AUG 


PARAM PARAM 
VY OUT<PARAM;B 
rf 


‘ENTER, IN ORDER, VALUES FOR THE FOLLOWING PARAMETERS: ' 


TRS 


"PANEL THICKNESS (INCHES) * 


B+ 


"RESIN COST (DOLLARS/POUND)' 


9 


B+B 0 


"WAX COST (DOLLARS/POUND)' 


"WOOD COST (DOLLARS/O0.D. POUND WAFERS)' 


r8 


B+B,0 


"FUEL COST (DOLLARS/MCF NATURAL GAS)' 


B+B 0 


ff 


‘OTHER VARIABLE COSTS (DOLLARS/MSF)' 


Srae: 


B+B8 0 


"SELLING PRICE AT THE FACTORY GATE (DOLLARS/MSF)' 


B+B 0 

Af 

OUT<B 

a 

I/O FUNCTION FOR SETTING UP USER PROVIDED 
INPUTS. 


ROW ROW 


ROW:(7~2+1 1,9w)pw 

A 

RESTRUCTURE w AS MATRIX WITH AT LEAST ONE ROW. 

(ONLY LAST TWO COORDINATES OF STRUCTURE ARE RETAINED.) 


>p>DvDD 
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xxxx  SEARCHMAX SEARCHMAX 


O| V FIRSTSTEP SEARCHMAX FIRSTBASE;OLDBASE;CURRENT ;STEPSAVE 

1| STATUS<0 7 
2| a STEPSAVE<STEPS 

Si) ZNITIALIZE: 

Wl a STEPS<STEPSAVE 

5| a MODE<O0 

6| BASEPOINT+FIRSTBASE 

IN STERSIZECRIRST STEP. 


8| CURRENT<FUNCTION BASEPOINT 


al@ | 9 
D1 | OI ORO I III III III I III I I IOI I A A A IC He 
’ 

ee M.C. RIGS DEN. P 
Biss 

AnSalig ete 

14| ‘STARTING POINT - PROCESS VARIABLES: *, 8 3 ¥(~14BASEPOINT ) 

15| ‘INITIAL VALUE OF LAGRANGE FUNCTION: ', 32 3 vFUNCTION BASEPOINT 

16| ‘LAMBDA IS: *, 32 3 ¥LAMBDA 

17| ‘VARIABLE REVENUE (DOLLARS/SHIFT): ', 32 3 ¥COST OBJFCN BASEPOTI 

104 MEE: BOND (PSI) IS: ‘, 32 3 ¥CONFCN ~14BASEPOINT 

19| START: 


20| TEST<CURRENT 

21| TEMPBASE+BASEPOINT 

22| +PRINTOUT IF(CONFCN ~14TEMPBASE)<MINIMUM 
23| EXPLOREMAX 

24| +PRINTOUT IF(CONFCN ~1+TEMPBASE)<MINIMUM 
25| +PATTERN IF TEST>CURRENT 

26| +PRINTOUT 


27| a*zPRINTOUT IF(+/STEPS(1;12STEPS(3;])=0 


OA 


tor oj Cl oO Aa ee , 


ae a - ee te 


ae 


; ee 


new 
5 oe 


Ve n : — : 
+ i] iors 
jt» 4) gta u gh ie ree » * 


* (heat ac teainny meta 


® oa 
7—F @&@ 


eo. isa, 
ante: ee 
vb cubs int 

vhowwdad be 

7, 

‘TTS. Vee 

' Ae. lef 
ene 
vanes 44 gf : 
oe | 
OOS exe 
ieee Dee 
“i ktees O56 
Te rien (OR 


. , Tid {ieee hay . 


28| ASTEPS( Ls leSTEPS (is IxSTEPSI 23) 

29| aSTEPSIZE<,STEPS(1;] 

30| aMODE+1 

31| axSTART 

32| PATTERN: 

33| +START IF(CONFCN ~14TEMPBASE.) <MINIMUM 

34| OLDBASE<BASEPOINT 

35| BASEPOINT+TEMPBASE 

36! CURRENT<TEST 

BFA) tat 

38 | ‘INTERMEDIATE - PROCESS VARIABLES: ", 8 3 WC 14BASEPOTINT ) 

39| ‘INTERMEDIATE VALUE OF LAGRANGE FUNCTION: ', 32 3 ¥FUNCTION BASEPOINT 
40| ‘VARIABLE REVENUE (DOLLARS/SHIFT): t, 32 3 ¥COST OBJFCN BASEPOL 
tay Dae BOND (PST) IS: 32) 3) ¥CONEGN, IV BASEPOINT. 
42| +SKIP IF STATUS=1 

43| TEMPBASE+(2xTEMPBASE)-OLDBASE 

44| TEST<FUNCTION TEMPBASE 

US SKIP: 

46| EXPLOREMAX 

47 | +PATTERN IF TEST>CURRENT 

48| +START 


49| PRINTOUT: 


SiO} au 

Si] Ae 

SD | 8 KI ORI I I II IIE II EI FO II I I I He 
kKkx«! 

53| ‘OPTIMAL SOLUTION FOR LAMBDA = ', 9 3 ¥LAMBDA 

Sales 

Gis} 9 Mee RaGe. DEN. Jee 
Je BASSS: 

56| 'FINAL POINT - PROCESS VARIABLES: ', 8 3 ¥(°1+BASEPOINT ) 


57| 'FINAL VALUE OF LAGRANGE FUNCTION: ', 32 3 ¥FUNCTION BASEPOINT 
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58| "VARIABLE REVENUE (DOLLARS/SHIFT): ', 32 3 ¥COST OBJFCN BASEPOI 
Sol ESE BOND (PSI) IS: ', 32 3 ¥CONFCN ~14BASEPOINT 
Goi} ov 
Balt 0% 


62| +OUT IFNOTO ‘WOULD YOU LIKE TO CONTINUE THIS ANALYSIS?! 
63| T+GAPQ 

64| +PRINTOUT IF T=1 

651 COSTQ 

66| BASEQ 

67| LAMBDAQ 1 

68| STATUSQ 

69| MINQ 1 


70)| SINITIALIZE 


Wak |p BG 

12) @OUnes 

73| tr 

TU] NK KK KK KKK KKK KK KKK END OF ANALYSIS KKK KK KKK KKK KKK KK OK 
xkkxt 

PS G 


76| a THE PRIMARY FUNCTION IN MAXPRESS. 

77| em SEARCHMAX PERFORMS THE CONSTRAINED OPTIMIZATION 
78| a OF THE WAFERBOARD PRESS CYCLE. 

79| a THIS FUNCTION SHOULD PROBABLY BE BROKEN DOWN 


80| a SO THAT I/0 COMPONENTS ARE SEPARATED. 
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SELECT SELECT 


V OUT<B SELECT A3;T3A131I3;COUNT 

BASE<B 

AtACAIAL31]3] 

A<A BY.iltpA 

T~(AL34]2MINIMUM)/11+ pA 

A<A[T3] 

A1l+(pA)p0 

I<2 

COUNT+1 

AVE RRS NE Seo | 

START: 

COUNT<COUNT+1 

>END IF COUNT>1+pA 

*START IF(AL1COUNT;3]>A(COUNT;3])(AL1COUNT 34J2A[COUNT;4)) 
A1[I3]+A(COUNT;] 

I<+I+1 

>START 

END: 

T+(A1[39]40)/riitpAl 

A1+A1(T3] 

A1+A1(¥A1(33]3] 

OUT<A1[13;] 

A 

THIS FUNCTION IS USED BY THE GAP SEARCH 
ROUTINE (GAPQ) TO SELECT THE ELEMENT 
WHICH PROVIDES THE GREATEST PAYOFF WHILE 
STILL MEETING PANEL QUALITY CONSTRAINTS. 


Dr DDD 


START START 


V START;A3;STEPSIZE;FIRSTBASE 


ts 
9) 9; 


‘RUN DATE *  yarsi ved 


bet) 


COST«<PARAM 

,f 

A+VARIABLES 

MIN<A[31] 

MAX+A[32] 

STEPSIZE+«A(33] 

‘STARTING POINT FOR VECTOR SEARCH;' 
FIRSTBASE<INITIAL,O 

LAMBDAQ@ 0 

MINQ 0 

STATUSQ 

STEPSET STEPSIZE 

STEPSIZE SEARCHMAX FIRSTBASE 
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STATUSQ 


VY STATUSQ 


tal | 


>SET IEY ESTO) "WOULD YOU “LIKE TO USE EXPEORATORY (MOV ES WONENG* 


STATUS<O 

>0 

SETH 

STATUS<1 

A 

AN I/O FUNCTION FOR CHANGING THE STATUS OF 
THE QUESTION POSED. 


BD D 


STEPSET 


WY SIMBA IIE RIP ASIOGDI 2 ISDA MSD IO SAGEM 
CONSTEP<0 

STEPS+(MAX-MIN)+STEPSIZE 

T<—yt 

SAM EES 

+OUT IF I>pSTEPS 
INT<MIN(I),MIN(IJ+STEPSIZE(IJxiSTEPS(TI]) 
CONSTEP<CONSTEP ON1 INT 

I<+I+1 

>START 

OUT: 

CONSTEP+ 1 0 +CONSTEP 

A 

THIS FUNCTION SETS UP THE STEP SIZE 
FOR THE CONSTRAINT DATA VARIABLE. 
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VARIABLES VARIABLES 


V OUT<VARIABLES:3B:3A3Q 

B< 4 3 p0 

e"¢ 

"ENTER, FOR EACH PROCESS VARIABLE, THE FOLLOWING PARAMETERS: ' 
"MINIMUM POSSIBLE VALUE, MAXIMUM POSSIBLE VALUE, STEP SIZE:' 
"xxxNOTE: STEP SIZE MUST EQUAL INTERVAL SIZE OF CONSTRAINT MATRIX.' 


as) 


‘EXAMPLE FOR MOISTURE CONTENT: ' 


Wy 


J i Eh oS 


ae 


"MOISTURE CONTENT OF THE PANEL (PERCENT)' 
ae 

AO] 

Q+11 

+(32%p,4)/ERROR 

B[1;1]<A 

°° 


‘RESIN CONTENT OF THE PANEL (PERCENT OF O.D. WEIGHT)' 
sf 

AQ) 

Q<18 

>(3%p,4)/ERROR 

BC2;]<A 

¢ 


‘PANEL DENSITY (PCF)' 
7a } 

A+D 

Q+25 

+(32p,4)/ERROR 
BL33;]<A 

26: 


‘PRESS TIME (MINUTES). 
IAS | 

A+] 

Q<+32 

>(34p,4)/ERROR 

BL4;]<A 

OUT<B 

>0 

ERROR: 


wae 


"WRONG ANSWER, PLEASE TRY AGAIN' 


9 


+ 
Le) 


AN I/O FUNCTION FOR ESTABLISHING VALUES FOR 
THE DECISION VARIABLES. 


RB DDD 
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13. APPENDIX III - GRAPHICAL REPRESENTATIONS OF CONSTRAINT 
DATA USED IN MAXPRESS 
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14. APPENDIX IV - SAMPLE RUN OF MAXPRESS 


NO: 


KOR RRR KKK KK ET RR RR KKK KK EE EK RK KEKE KT RE KK 


KR RK KKK KK KK RRR KK RR KK REE KEK KEKE EET EEK 


**  WAFERBOARD PRODUCTION OPTIMIZATION MODEL 


se ALBERTA RESEARCH COUNCIL 

is FOREST PRODUCTS PROGRAM: FP-19 
ee AUTHOR: TOM I. GRABOWSKI 

ars LAST REVISION: APRIL 29, 1982 


*** 
** 
KO 
ok 


** 


OR RK OK KOK KK KKK KO KKK RR ERK KEK KKK RK KKK KEK KKK 


KKK KK KKK KKK KKK KKK KR KKK ERK KK KEK KKK KKK RK EK KEK KEK 


*****WARNING: YOUR DATA FILE MUST BE COPIED INTO THIS WORKSPACE . 


ENTER ONE OF THE FOLLOWING KEYWORDS: 

EXPLAIN (GIVES A DESCRIPTION OF THE PROGRAM) 
START (INITIATE RUNNING OF THE PROGRAM) 
START 


RUN DATE TIS 25 SaaS S 


ENTER, IN ORDER, VALUES FOR THE FOLLOWING PARAMETERS: 


PANEL THICKNESS (INCHES) 
.4375 

RESIN COST (DOLLARS/POUND ) 
17. 

WAX COST (DOLLARS/POUND ) 
BOw 

WOOD COST (DOLLARS/O.D. POUND WAFERS) 
BOS 

FUEL COST (DOLLARS/MCF NATURAL GAS) 
15 

OTHER VARIABLE COSTS (DOLLARS/MSF ) 


50 


SELLING PRICE AT THE FACTORY GATE (DOLLARS/MSF ) 


200 


108 


(cae =f n 

e ae — » 
peiareia 

ye >. pad ral ' 


s 


rep eae e445 ’ 


vito" 93 2am ween 


aw 3 WY clea a a r 
thd pre 
ar? 4¢ TAS, 


—~ 
‘ 
g + ASE 
we _ 
eels x = 7 
a "(1 8) 
by 


soa Llane aa at 


7 
te: oy 
24,000) Teenie a 
: a a 
4500) oe 
ery 
’ 7 : if 7 > > 7 Taipew 
q i“ 


pine | (4: 7. uo #6 42068 ; 


104 


ENTER, FOR EACH PROCESS VARIABLE, THE FOLLOWING PARAMETERS: 
MINIMUM POSSIBLE VALUE, MAXIMUM POSSIBLE VALUE, STEP SIZE: 
***NOTE: STEP SIZE MUST EQUAL INTERVAL SIZE OF CONSTRAINT MATRIX. 
EXAMPLE FOR MOISTURE CONTENT: 
@ & .& 
MOISTURE CONTENT OF THE PANEL (PERCENT) 
Ge) 2S 
RESIN CONTENT OF THE PANEL (PERCENT OF O.D. WEIGHT) 
AL ee shet Gy VERS 
PANEL DENSITY (PCF) 
32 48 2 
PRESS TIME (MINUTES) 
Si Jiao 
STARTING POINT FOR VECTOR SEARCH; 
ENTER, IN ORDER, VALUES FOR EACH PROCESS VARIABLE: 
MOISTURE CONTENT (PERCENT) 
7 
RESIN CONTENT (PERCENT) 
22, PANO) 
PANEL DENSITY (PCF) 
40 
PRESS TIME (MINUTES) 


5 


ENTER THE VALUE FOR LAMBDA 


150 


ENTER VALUE FOR MINIMUM INTERNAL BOND 


O 


WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 
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REC 


STARTING POINT - PROCESS VARIABLES: 
INITIAL VALUE OF LAGRANGE FUNCTION: 
LAMBDA IS: 

VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (ODOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


7 


M_G. 


. 000 


. 000 


. 500 


. 500 


2b. 


250 


SIRS) 


Ye 


» BAO 


DEN. 


40.000 5 


14025. 
1SOF 
5540. 
56. 


38 .000 4. 
14758. 
6878. 

eA 


36 . 000 4 
15071 


8910. 
.068 


41 


34.000 4 


32 


PRESS 


. 000 


556 
000 
360 
568 


500 
844 
821 
533 


. 000 
RON 


858 


. 000 
15086. 
10150. 
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026 
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OPTIMAL SOLUTION FOR LAMBDA = 150.000 

MC. RECZ DEN. PRESS 
FINAL POINT - PROCESS VARIABLES: TOO 2.250 34.000 4.000 
FINAL VALUE OF LAGRANGE FUNCTION: 15086 .437 
VARIABLE PROFIT (DOLLARS/SHIFT): 10150.026 
INTERNAL BOND (PSI) IS: 32.909 
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WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 
WOULD YOU LIKE TO PERF.ORM A GAP SEARCH? N 
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WOULD YOU LIKE TO CHANGE THE VARIABLE COST FIGURES? N 
WOULD YOU LIKE TO CHANGE THE STARTING POINT FOR THE SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VALUE OF LAMBDA? Y 


ENTER THE VALUE FOR LAMBDA 


160 


WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 
WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND? N 


RK KK ROK ROKR KO RK KOK K KK KOK KE KR ROKK KR ROR OK RR KK RR KKK KK EK KEK KEK KKK KKK KK KOK KK 


STARTING POINT - PROCESS VARIABLES: 
INITIAL VALUE OF LAGRANGE FUNCTION: 
LAMBDA IS: 

VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


M.C. 


7.000 


7.000 


[500 


7.000 


R.C. DEN. PRESS 


26 250m 402 O00 5.000 
W499 18 236 
160.000 
5540.360 
56.568 

Ph B\ifia) ASKS} KOLOK0} 4.500 
15284.179 

6878 .821 

S}P4 ASSIS! 

ASUS) BS. CO, 4.500 
WiSreXs) , AAS 
7920.762 
47.547 

22 Ss} SG}. (CLON@) 4.500 
SSIS HS}, SISI7/ 
7810.251 
48.282 
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OPTIMAL SOLUTION FOR LAMBDA = 160.000 


FINAL POINT - PROCESS VARIABLES: 
FINAL VALUE OF LAGRANGE FUNCTION: 
VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


ReCEr DEN. PRESS 
Ph Sila) SKE} OO) 4.500 
SSS 5m Si 

TEES KO) XS) t, 

48.282 
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WOULD 
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CONTINUE THIS ANALYSIS? Y 
PERFORM A GAP SEARCH? N 


CHANGE THE VARIABLE COST FIGURES? N 
CHANGE THE STARTING POINT FOR THE SEARCH? N 


CHANGE THE VALUE OF LAMBDA? Y 
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ENTER THE VALUE FOR LAMBDA 


155 


WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 
WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND? N 
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STARTING POINT - PROCESS VARIABLES: 
INITIAL VALUE OF LAGRANGE FUNCTION: 
LAMBDA IS: 

VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (ODOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


M.C. Ra Ce DEN. PRESS 


7.000 2.250 40.000 5.000 
14308 .396 
155.000 
5540.360 
56.568 

7.000 25370 36 OOO 4.500 
150212 Sil2 

6878 .821 
22-993 

7.500 2 SUS SKS) OKO) 4.500 
15290.493 
7920.762 
47.547 

7.000 253 OSCE OOO 4.500 
15293 .947 
SMOR2on 


48.282 
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OPTIMAL SOLUTION FOR LAMBDA = 155.000 


FINAL POINT - PROCESS VARIABLES: 
FINAL VALUE OF LAGRANGE FUNCTION: 
VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


MGr RinGu DEN. PRESS 
7.000 2 SUS) SKS} (COXS) 4.500 
15293 .947 
TEVOR2 oN 


48.282 
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WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 

WOULD YOU LIKE TO PERFORM A GAP SEARCH? N 

WOULD YOU LIKE TO CHANGE THE VARIABLE COST FIGURES? N 

WOULD YOU LIKE TO CHANGE THE STARTING POINT FOR THE SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VALUE OF LAMBDA? Y 

ENTER THE VALUE FOR LAMBDA 


Syd 6S} 


WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 
WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND? N 
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M.C. IC. DEN. PRESS 
STARTING POINT - PROCESS VARIABLES: 7.000 2.250 40.000 5.000 
INITIAL VALUE OF LAGRANGE FUNCTION: 14166.976 
LAMBDA IS: Ksi22 . SKOXO) 
VARIABLE PROFIT (DOLLARS/SHIFT): 5540.360 
INTERNAL BOND (PSI) IS: 56.568 
INTERMEDIATE - PROCESS VARIABLES: 7 .OQO0O0 231 Io 1 OOO 4.500 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 14890.178 
VARIABLE PROFIT (DOLLARS/SHIFT): 6878. 821 
INTERNAL BOND (PSI) IS: S2e1919 3 
INTERMEDIATE - PROCESS VARIABLES: 72500 PL BIS) ENS. CKO) 4.000 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 15173 .686 
VARIABLE PROFIT (DOLLARS/SHIFT): 8910.858 
INTERNAL BOND (PSI) IS: 41.068 
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OPTIMAL SOLUTION FOR LAMBDA = 1522500 

MENGE R.C. DEN. ~ TPARIE SSS 
FINAL POINT - PROCESS VARIABLES: 25 OO Ph Sify B15), OO 4.000 
FINAL VALUE OF LAGRANGE FUNCTION: Wel 7s) GES 
VARIABLE PROFIT (DOLLARS/SHIFT): 8910.858 
INTERNAL BOND (PSI) IS: 41.068 


WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 
WOULD YOU LIKE TO PERFORM A GAP SEARCH? Y 
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KKK KEK KKK EK KKK KK Oe GAP SEARCH ROUTINE 
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ENTER, IN ORDER, VALUES FOR EACH PROCESS VARIABLE: 


MOISTURE CONTENT (PERCENT) 
Pos 

RESIN CONTENT (PERCENT) 
An iS) 

PANEL DENSITY (PCF) 
36 

PRESS TIME (MINUTES) 


4 


ENTER VALUE FOR MINIMUM INTERNAL BOND 


42 


ENTER THE VALUE FOR LAMBDA 


O22 


KEKE KKK K KK 


ENTER THE VALUE FOR PERTURBATION DEPTH (MAXIMUM IS THREE) 


1 
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OPTIMAL SOLUTION FOR LAMBDA = 


FINAL POINT - PROCESS VARIABLES: 


FINAL VALUE OF LAGRANGE FUNCTION: 
VARIABLE PROFIT (DOLLARS/SHIFT): 


INTERNAL BOND (PSI) IS: 


MEGr 
7 .O0O0 


RiGy DEN. PRESS 
2 DOOM s6r OOO 4.000 
2083.237 

8576 .022 

42.576 
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WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 
WOULD YOU LIKE TO PERFORM A GAP SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VARIABLE COST FIGURES? Y 
ENTER, IN ORDER, VALUES FOR THE FOLLOWING PARAMETERS: 
PANEL THICKNESS (INCHES) 
.4375 
RESIN COST (DOLLARS/POUND ) 
1.00 
WAX COST (DOLLARS/POUND ) 
oy/ 
wOOD COST (DOLLARS/O.D. POUND WAFERS) 
.O5 
FUEL COST (DOLLARS/MCF NATURAL GAS) 
1.5 
OTHER VARIABLE COSTS (DOLLARS/MSF ) 
50 


SELLING PRICE AT THE FACTORY GATE (DOLLARS/MSF ) 


200 


WOULD YOU LIKE TO CHANGE THE STARTING POINT FOR THE SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VALUE OF LAMBDA? N 

WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 

WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND? Y 

ENTER VALUE FOR MINIMUM INTERNAL BOND 
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M.C. Reice DEN. PRESS 
STARTING POINT - PROCESS VARIABLES: 7.000 2.250 40.000 5.000 
INITIAL VALUE OF LAGRANGE FUNCTION: 1259 8m 552 
LAMBDA IS: 152.500 
VARIABLE PROFIT (DOLLARS/SHIFT): SE)" i}. SSIS 
INTERNAL BOND (PSI) IS: 56.568 
INTERMEDIATE - PROCESS VARIABLES: 7.000 AX SIfis) Ss} (O10) 5.000 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: USS o 1°72 
VARIABLE PROFIT (DOLLARS/SHIFT): 4560. 266 
INTERNAL BOND (PSI) IS: 56.338 
INTERMEDIATE - PROCESS VARIABLES: 7.000 2S OS On OOO 4.500 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 13456.744 
VARIABLE PROFIT (DOLLARS/SHIFT): GOSSruS 
INTERNAL BOND (PSI) IS: 48.282 
INTERMEDIATE - PROCESS VARIABLES: 7.500 2H DOMES OOOO 4.500 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 13486 .484 
VARIABLE PROFIT (DOLLARS/SHIFT): 6487 .967 
INTERNAL BOND (PSI) IS: 45.892 
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OPTIMAL SOLUTION FOR LAMBDA = 152.500 

MECe REC. DEN. PRESS 
FINAL POINT - PROCESS VARIABLES: 7.500 22250) 736 000 4.500 
FINAL VALUE OF LAGRANGE FUNCTION: 13486 .484 
VARIABLE PROFIT (DOLLARS/SHIFT): 6487 .967 
INTERNAL BOND (PSI) IS: 45.892 


WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 
WOULD YOU LIKE TO PERFORM A GAP SEARCH? Y 
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ENTER, IN ORDER, VALUES FOR EACH PROCESS VARIABLE: 
MOISTURE CONTENT (PERCENT ) 
mais) 
RESIN CONTENT (PERCENT ) 
229 
PANELS DENS TTY. (PCE) 
36 


PRESS TIME (MINUTES) 


4.5 


ENTER VALUE FOR MINIMUM INTERNAL BOND 


42 


ENTER THE VALUE FOR LAMBDA 


WS oS) 


ENTER THE VALUE FOR PERTURBATION DEPTH (MAXIMUM IS 3) 


3 


ROR Ke ROR ORR KK KKK KK RK KR KK EKKO KKK KEK KR KK RE KKK KR KEK KKK KR KKK KK KK K 


OPTIMAL SOLUTION FOR LAMBDA = SAD. KOKO) 

M.C. RC. DEN. PRESS 
FINAL POINT - PROCESS VARIABLES: tie SOO 2a OO OLOOO 4.500 
FINAL VALUE OF LAGRANGE FUNCTION: 295.674 
VARIABLE PROFIT (DOLLARS/SHIFT): 6762.934 


INTERNAL BOND (PSI) IS: 42.408 
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WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 
WOULD YOU LIKE TO PERFORM A GAP SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VARIABLE COST FIGURES? Y 
ENTER, IN ORDER, VALUES FOR THE FOLLOWING PARAMETERS: 
PANEL THICKNESS (INCHES) 
.4375 
RESIN COST (DOLLARS/POUND ) 
120 
WAX COST (DOLLARS/POUND ) 
FOW: 
WOOD COST (DOLLARS/O.D. POUND WAFERS) 
.05 
FUEL COST (DOLLARS/MCF NATURAL GAS) 
1.5 
OTHER VARIABLE COSTS (DOLLARS/MSF ) 
25 


SELLING PPICE AT THE FACTORY GATE (DOLLARS/MSF ) 


150 


WOULD YOU LIKE TO CHANGE THE STARTING POINT FOR THE SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VALUE OF LAMBDA? Y 
ENTER THE VALUE FOR LAMBDA 


WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 
WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND? Y 
ENTER VALUE FOR MINIMUM INTERNAL BOND 
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M.C. Rw. DEN. PRESS 
STARTING POINT - PROCESS VARIABLES: 7.000 2.250 40.000 5.000 
INITIAL VALUE OF LAGRANGE FUNCTION: 5316.654 
LAMBDA IS: 90.000 
VARIABLE PROFIT (DOLLARS/SHIFT): 225.536 
INTERNAL BOND (PSI) IS: 56.568 
INTERMEDIATE - PROCESS VARIABLES: 7.500 2.250 38.000 5.500 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 6180.475 
VARIABLE PROFIT (DOLLARS/SHIFT): 1134.881 
INTERNAL BOND (PSI) IS: 56.062 
INTERMEDIATE - PROCESS VARIABLES: 7.500 2.250 36.000 6.500 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 6659.785 
VARIABLE PROFIT (DOLLARS/SHIFT): 1655.977 
INTERNAL BOND (PSI) IS: 55.598 
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OPTIMAL SOLUTION FOR LAMBDA = 90.000 

Mace RaGr DEN. PRESS 
FINAL POINT - PROCESS VARIABLES: 7.500 2.250 36.000 6.500 
FINAL VALUE OF LAGRANGE FUNCTION: 6659.785 
VARIABLE PROFIT (DOLLARS/SHIFT): NODS) SAT 
INTERNAL BOND (PSI) IS: hs). SISK} 


WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 

WOULD YOU LIKE TO PERFORM A GAP SEARCH? N 

WOULD YOU LIKE TO CHANGE THE VARIABLE COST FIGURES? N 

WOULD YOU LIKE TO CHANGE THE STARTING POINT FOR THE SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VALUE OF LAMBDA? Y 

ENTER THE VALUE FOR LAMBDA 
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WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 
WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND? N 


yi 


7 °i pe Pe 
=| 
ot @ 


7 = ale 
oer 

’ a : 

ce er io. 


= ¥ ag cei 
=> Vee 
os } 
“56 _ ‘ps ahaa _ it 


va’) 
= a] FoF) we Pe’ v1 7 
ts cT)* See iar 


) OF) grove rg 


4 
- - 
-_ 


ee ee | 
ear Ye 
-? @ae 


= a “SS 9 iC 
_|° ornare ams vl 
> son Bey x 


UL ead 


is ov Guan 


1ahe ioG@) ; 
a. 


TiS 


RO OR OR CR RO ROR OK OR RK KK ROKR RT RKTT KR RK KET EEK EEK KRE KKK KK KK EE 


M.C. 


STARTING POINT - PROCESS VARIABLES: 
INITIAL VALUE OF LAGRANGE FUNCTION: 
LAMBDA IS: 

VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE GF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 


INTERMEDIATE - PROCESS VARIABLES: 


INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 


VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 
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PRESS 


5.000 


5146 .950 
87.000 
225.536 
56.568 


5.500 


6012 .288 
1134.881 
DomOo2 


5.500 


655m o2 
2779.247 
43.356 


5.000 


6647 .408 
SOSiemin 2: 
41.267 


4.000 


7107 .088 
5226 .934 
Zonal 


3.000 


SO). OV 
7336 .987 


1.940 


3.000 


7646 .824 
7490.224 


1.800 
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OPTIMAL SOLUTION FOR LAMBDA = 87.000 


FINAL POINT - PROCESS VARIABLES: 
FINAL VALUE OF LAGRANGE FUNCTION: 
VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BONO (PSI) IS: 
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PERFORM A GAP SEARCH? 
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CONTINUE THIS ANALYSIS? 
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CHANGE THE VARIABLE COST FIGURES? N 
CHANGE THE STARTING POINT FOR THE SEARCH? N 
CHANGE THE VALUE OF LAMBDA? Y 


LIKE TO USE EXPLORATORY MOVES ONLY? N 
LIKE TO CHANGE MINIMUM INTERNAL BOND? N 
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STARTING POINT - PROCESS VARIABLES: 
INITIAL VALUE OF LAGRANGE FUNCTION: 
LAMBDA IS: 

VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 
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OPTIMAL SOLUTION FOR LAMBDA = 88.000 


FINAL POINT - PROCESS VARIABLES: 
FINAL VALUE OF LAGRANGE FUNCTION: 
VARIABLE PROFIT (DOLLARS/SHIFT): 
INTERNAL BOND (PSI) IS: 
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WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 

WOULD YOU LIKE TO PERFORM A GAP SEARCH? N 

WOULD YOU LIKE TO CHANGE THE VARIABLE COST FIGURES? N 

WOULD YOU LIKE TO CHANGE THE STARTING POINT FOR THE SEARCH? N 
WOULD YOU LIKE TO CHANGE THE VALUE OF LAMBDA? Y 

ENTER THE VALUE FOR LAMBDA 
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WOULD YOU LIKE TO USE EXPLORATORY MOVES ONLY? N 
WOULD YOU LIKE TO CHANGE MINIMUM INTERNAL BOND? Y 


ENTER VALUE FOR MINIMUM INTERNAL BOND 
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M.C. R.C. DEN. PRESS 
STARTING POINT - PROCESS VARIABLES: 7.000 2.250 40.000 5.000 
INITIAL VALUE OF LAGRANGE FUNCTION: 5146.950 
LAMBDA IS: 87.000 
VARIABLE PROFIT (DOLLARS/SHIFT): 225.536 
INTERNAL BOND (PSI) IS: 56.568 
INTERMEDIATE - PROCESS VARIABLES: Uf 5 SOO) 2.250 38.000 5.500 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 6012.288 
VARIABLE PROFIT (DOLLARS/SHIFT): 1134.881 
INTERNAL BOND (PSI) IS: 56.062 
INTERMEDIATE - PROCESS VARIABLES: 7.500 2.250 34.000 B) 5 BIOXe) 
INTERMEDIATE VALUE OF LAGRANGE FUNCTION: 6551.262 
VARIABLE PROFIT (OOLLARS/SHIFT): PRT TE) PY 
INTERNAL BOND (PSI) IS: 43.356 
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OPTIMAL SOLUTION FOR LAMBDA = 87.000 

MEGr. RG DEN. PRESS 
FINAL POINT - PROCESS VARIABLES: Tg 5 OO 2.250 34.000 5.500 
FINAL VALUE OF LAGRANGE FUNCTION: 6551.262 
VARIABLE REVENUE (DOLLARS/SHIFT): 2779 .247 
INTERNAL BOND (PSI) IS: 43.356 


WOULD YOU LIKE TO CONTINUE THIS ANALYSIS? Y 
WOULD YOU LIKE TO PERFORM A GAP SEARCH? Y 
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ENTER, IN ORDER, VALUES FOR EACH PROCESS VARIABLE: 
MOISTURE CONTENT (PERCENT) 

Uf 332) 
RESIN CONTENT (PERCENT ) 

225 
PANEL DENSITY (PCF) 

36 
PRESS TIME (MINUTES) 


625 


ENTER VALUE FOR MINIMUM INTERNAL BOND 
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ENTER THE VALUE FOR LAMBDA 
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ENTER THE VALUE FOR PERTURBATION DEPTH (MAXIMUM IS THREE) 
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OPTIMAL SOLUTION FOR LAMBDA = 87.000 

M.C. RECs DEN. PRESS 
FINAL POINT - PROCESS VARIABLES: 6.500 2.250 34.000 5.000 
FINAL VALUE OF LAGRANGE FUNCTION: 6519.939 
VARIABLE PROFIT (DOLLARS/SHIFT): 2855 .678 
INTERNAL BOND (PSI) IS: ADRS 
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